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SUMMARY 


An  experimental  and  analytical  program  was  conducted  to  develop  a 
methodology  to  predict  total  recession,  including  the  contribution  by 
particulate  impact  erosion,  in  solid  propellant  rocket  nozzles.  The  program 
was  divided  into  two  phases.  Under  Phase  1,  a  methodology  was  developed  for 
predicting  recession  due  to  particle  impact  erosion  in  the  expansion  (exit 
cone)  region  of  nozzles.  This  effort  was  performed  by  the  Aerojet  Strategic 
Propulsion  Company  (ASPC)  under  a  subcontract  to  the  Aerotherm  Division  of 
Acurex  Corporation.  Under  Phase  II,  a  methodology  was  developed  for 
predicting  recession  due  to  erosion  coupled  with  ablation  in  the  entrance  and 
throat  regions  of  nozzles.  This  effort  was  performed  by  Aerotherm. 

Results  from  this  program  are  reported  in  four  separate  volumes,  as 
follows; 


Volume  I: 


Volume  II: 


Volume  III: 


Volume  IV; 


Recession  Prediction  Methodology  for  Rocket  Nozzle 
Entrance  and  Throat  Regions  (Reference  1) 

Erosion  Model  Development  (Expansion  Region) 
(Reference  2) 

User's  Manual  —  Chamber  Flowfield  Code  (CFC)  and 
Charring  Material  Ablation-Erosion  (CMAE)  Code 
User's  Manual  —  Erosion  Prediction  Procedure  for 
Rocket  Nozzle  Expansion  Region  (Reference  3) 
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The  remainder  of  this  summary  provides  a  brief  overview  of  the 
recession  prediction  methodologies  developed  for  the  nozzle  expansion  region 


(reported  in  Volume  II)  and  for  the  nozzle  entrance  and  throat  regions 
(reported  in  Volume  I). 

Nozzle  Expansion  Region 

A  procedure  was  developed  for  predicting  the  erosion  that  occurs  in  a 
solid  rocket  nozzle  exit  cone  when  aluminum  oxide  (AI2O3)  particles  in  the 
flow  impact  the  nozzle  wall  in  the  supersonic  region  downstream  of  the  nozzle 
throat.  This  procedure  consists  of  a  two-phase  flow  solution  and  two 
independent  mechanical  erosion  models.  These  models  calculate  the  nozzle  wall 
erosion  patterns  in  carbon-carbon  and  bulk  graphite  materials  as  a  function  of 
the  impinging  particle  flow  properties,  the  impact  angle,  and  the  nozzle  wall 
material  properties. 

An  empirical  "G-law"  erosion  model  was  developed  for  nozzle  expansion 
regions  from  Dust  Erosion  Tunnel  (DET)  test  data  and  small  solid  motor  firing 
data.  Specimens  of  carbon-carbon,  bulk  graphite,  and  tantalum/tungsten  were 
exposed  to  particle-laden  test  gases.  A  debris  layer  shielding  erosion  model 
was  also  developed  for  the  expansion  region.  In  contrast  to  the  empirical 
G-law,  this  model  was  formulated  to  represent  the  fundamental  physics  of  the 
erosion  process. 

The  TD2P  two-phase  flow  code  contained  in  the  Interim  Solid  Performance 
Program  (ISPP)  was  modified  to  incorporate  the  above  two  erosion  models.  This 
version  of  T02P  constitutes  the  methodology  for  predicting  particle  impact 
conditions  at  the  walls  of  nozzle  exit  cones  and  the  resulting  wall  erosion. 
Predictions  were  made  using  both  erosion  models  for  various  solid  rocket  motor 
test  firings  where  erosion  occurred.  Predictions  using  the  G-law  model  were 
generally  acceptable  relative  to  test  results,  providing  the  particle  flow  and 


6 


impact  conditions  were  within  the  range  covered  by  the  DET  and  small  motor 
experiments.  Predictions  using  the  debris  layer  model  were  also  generally 
satisfactory,  but  not  as  accurate  as  G-law  predictions  when  particle  impact 
conditions  were  within  the  DET  and  small  motor  test  range.  However,  the 
debris  layer  model  was  the  most  accurate  for  conditions  outside  the 
experimental  range. 

Nozzle  Entrance  and  Throat  Regions 

A  methodology  was  developed  for  predicting  the  total  recession, 
including  particle  impact  erosion,  for  the  entrance  and  throat  regions  of 
solid  rocket  nozzles.  Major  program  efforts  included  the  development  of  a 
motor  chamber-nozzle  flowfield  code  and  a  model  for  describing  the  erosion 
process  in  the  subsonic  regions  of  the  nozzle. 

The  Fully  Coupled  Transonic  (FCT)  two-phase  flow  module  of  the  Solid 
Performance  Prediction  (SPP)  program  (Reference  4)  was  modified  to  compute 
particle  trajectory  crossovers  and  particle  impact  conditions  at  the  nozzle 
wall.  Also,  the  calculation  domain  was  extended  forward  to  account  for  the 
motor  chamber  and  propellant  grain,  including  mass  injection  from  the  grain. 
However,  flow  within  the  submerged  region  of  nozzles  remains  to  be  incluued. 
This  modification  is  referred  to  as  the  Chamber  Flowfield  Code  (CFC). 

Empirical  subsonic  erosion  models  were  developed  for  predicting 
recession  due  to  the  impact  of  aluminum  (Al)  and  AI2O3  particles.  Models  were 
constructed  in  the  form  of  a  G-law  with  two  components:  mechanical  erosion 
due  to  mechanical  mass  removal  and  chemical  erosion  due  to  particle-surface 
chemical  reactions.  However,  the  models  are  based  on  bulk  graphite  data  and 
provide  only  rough  estimates  for  the  erosion  of  carbon-carbon  composites.  The 
Charring  Material  Ablation  (CMA)  code  (Reference  5)  was  modified  to  include 
t^ese  model'  in  predictions  of  recession  resulting  from  ablation  and  erosion. 
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This  modification  is  designated  the  Charring  Material  Ablation-Erosion  (CMAE) 
code. 

The  recession  prediction  methodology  for  subsonic  flow  regimes  wa^ 
demonstrated  for  two  motor  firings  having  suspected  particle  erosion  in  the 
entrance  region  of  carbon-carbon  composite  nozzles.  The  demonstration  results 
substantiated  the  prediction  of  the  onset  of  particle  impact  for  the  two 
motors.  However,  recession  predictions  using  the  bulk  graphite  erosion  model 
provided  only  a  rough  estimate  of  material  removed. 

Although  not  currently  included  in  the  recession  prediction 
methodology,  particle  breakup  in  nozzle  entrance  and  throat  regions  was 
modeled  in  terms  of  Bond  number  and  nondimensional  time.  This  model  was 
based  upon  cold- flow  data,  and  its  applicability  to  hot- firing  conditions  has 
not  been  demonstrated. 
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SECTION  1 


INTRODUCTION 


Described  herein  are  the  input  and  output  features  of  two  analytical 
procedures  developed  under  the  Particle  Impact  Erosion  Program:  The  Chamber 
Flowfield  Code  (CFC)  and  the  Charring  Material  Ablation-Erosion  (CMAE)  code. 
These  two  codes  are  key  constituents  of  the  methodology  for  predicting  total 
recession,  including  the  effects  of  particulate  erosion,  in  the  subsonic  and 
transonic  regions  of  solid  propellant  rocket  nozzles.  Details  of  the 
methodology  are  given  in  Volume  I  of  the  Particle  Impact  Erosion  Program  Final 
Report  (Reference  1).  This  report  is  divided  into  two  parts,  as  follows. 

1.1  PART  1  —  CHAMBER  FLOWFIELD  CODE 

Part  1  is  a  user's  manual  for  the  Chamber  Flowfield  Code  (CFC),  which 
calculates  particle  impact  conditions  in  the  entrance  and  throat  regions  of 
solid  propellant  rocket  nozzles.  CFC  is  an  extension  of  the  Fully  Coupled 
Transonic  (FCT)  module  of  the  Solid  Performance  Program  (SPP).  It  is  presumed 
that  the  user  is  thoroughly  familiar  with  the  use  of  the  FCT  module  of  the  SPP 
code.  Since  the  basic  features  of  the  FCT  module  are  preserved  in  CFC,  this 
manual  should  be  used  in  conjunction  with  the  SPP  user's  manual  (Reference  4). 
Only  modifications  to  the  SPP  input-output  procedure  are  detailed  herein. 

Major  modifications  made  to  the  FCT  code  as  part  of  the  CFC  development  are 


described  below. 


Geometry 

The  calculation  domain  was  extended  forward  into  the  chamber  to  account 
for  the  effects  of  chamber  and  grain  configuration  on  particle  impact 
conditions.  As  in  FCT,  only  axisymmetric  geometries  are  allowed,  and  the  flow 
within  the  submerged  region  near  the  backside  of  the  nozzle  is  not  computed. 
However,  the  effects  of  the  flow  from  the  submerged  region  are  treated  in  an 
approximate  manner.  Also,  complex  prepellant  geometries,  such  as  those  with 
star  grains  and  slots,  can  only  be  treated  approximately  with  the  code. 
Boundary  Conditions 

The  boundary  conditions  were  modified  to  handle  a  specified  mass  flux 
and/or  a  grain  burn-rate  distribution  along  the  input  chamber-grain  boundary. 
If  the  chamber  has  a  front-end  closure  (with  no  mass  injection),  the  code  can 
treat  that  boundary  separately  using  slip-flow  conditions. 

Particle  Flow  Calculations 

The  particle  trajectory  calculations  were  modified  to  allow  for  impact 
at  the  wall,  particle  trajectory  crossover,  and  particles  crossing  the  nozzle 
axis.  Locations  of  particle  impact  on  the  wall  are  identified  for  each 
particle  size  group,  and  the  distributions  of  particle  impingement  velocity, 
angle,  density  and  mass  flux  along  the  wall  are  output.  The  particle  size 
distribution  calculations  were  also  modified.  For  erosion  calculations  the 
large  aluminum  particles  coming  off  the  grain  are  of  more  interest  than  the 
smaller  aluminum  oxide  particles  formed  in  the  chamber.  Hence,  CFG  includes 
an  option  to  calculate  the  mean  diameter  and  size  distribution  of  the  aluminum 
agglomerates  leaving  the  grain  surface.  Presently,  as  in  FCT,  the  CFG  code 
does  not  account  for  particle  size  change  mechanisms. 
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1.2  PART  2  —  CHARRING  MATERIAL  ABLATION-EROSION  (CMAE)  CODE 

Part  2  is  a  user's  manual  for  the  Charring  Material  Ablation-Erosion 
(CMAE)  code.  This  code  calculates  the  total  recession,  including  erosion  due 
to  particle  impact,  and  one-dimensional  in-depth  transient  conduction  in 
rocket  nozzle  materials.  CMAE  is  an  extension  of  the  Charring  Material 
Ablation  (CMA)  code  and  preserves  the  basic  features  of  CMA.  Only  the 
modifications  to  the  CMA  input-output  procedure  are  detailed  here;  therefore, 
this  manual  should  be  used  in  conjunction  with  the  CMA  user’s  manual 
(Reference  5). 

The  major  modification  to  CMA  is  the  inclusion  of  erosion  in  the 
recession  calculations.  CMAE  contains  an  erosion  model,  developed  in  the 
Particle  Impact  Erosion  program,  for  the  subsonic-transonic  regions  of  rocket 
nozzles.  Details  of  the  model  development  are  given  in  Reference  1.  In 
summary,  the  erosion  model  contains  two  components:  mechanical  erosion  due  to 
mechanical  mass  removal  of  material,  and  chemical  erosion  due  to  chemical 
reactions  between  molten  particles  and  the  surface.  Each  erosion  component  is 
in  the  form  of  a  G-law,  G  being  the  ratio  of  mass  flux  of  target  material 
removed  to  mass  flux  of  particles  incident  upon  the  surface.  G  is  related  to 
particle  impact  conditions  in  the  form  of  power  laws  commonly  used  in  reentry 
vehicle  erosion  research  (Reference  11). 

Following  sections  briefly  describe  the  CFC  and  CMAE  codes  and  the 
changes  to  the  SPP  and  CMA  user's  manuals  required  to  use  them.  Sample  cases 
executed  with  CFC  and  CMAE  are  also  presented. 
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PART  1 

CHAMBER  FLOWFIELD  CODE  (CFC) 
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SECTION  2 


CFC  ANALYTICAL  BASIS 

Particle  impact  conditions  can  be  determined  if  the  particle 
trajectories  within  the  nozzle  are  known.  Since  the  particle  and  gas  phases 
affect  one  another,  a  coupled  two-phase  flowfield  calculation  must  be  made  to 
establish  the  particle  trajectories.  The  coupled  calculation  must  also 
account  for  other  important  nozzle  flow  characteristics,  such  as 
compressibility,  two-dimensional  (axisymmetric)  flow,  and  the  presence  of 
different  sized  particles  within  the  nozzle.  Additionally,  since  the  particle 
trajectories  in  the  entrance  region  are  expected  to  be  strongly  influenced  by 
the  chamber  and  grain  geometry,  the  flow  calculation  must  include  most,  if  not 
all,  of  the  chamber  region,  as  well  as  the  nozzle  entrance  and  throat  regions. 
To  achieve  these  requirements  in  the  i/C,  three  major  changes  were  made  to  the 
FCT  code: 

•  Extension  of  the  calculation  domain  to  include  the  chamber 

•  Modification  of  inflow  boundary  conditions  to  account  for  mass 
injection  from  the  propellant  grain 

•  Modification  to  particle  trajectory  calculations  to  evaluate 
trajectory  crossovers  and  particle  impact  at  the  wall 

All  changes  were  ma,de  to  existing  subroutines;  no  new  subroutines  were 
added.  The  changes  are  discussed  in  the  subsection  below,  where  it  is 
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assumed  that  the  reader  is  well  acquainted  with  the  FCT  module  of  SPP  and  its 
documentation  (Reference  4). 

2.1  GEOMETRY 

The  calculation  domain  of  the  base  (FCT)  code  includes  only  the  nozzle 
entrance  and  throat  regions.  The  starting  line  for  the  calculations  is  a 
fictitious  inflow  boundary  drawn  in  the  vicinity  of  the  nozzle  inlet.  Since 
the  effect  of  the  chamber  and  grain  geometry  is  expected  to  be  significant  on 
particle  impact  conditions,  the  calculation  domain  was  extended  into  the 
chamber. 

The  FCT  code  uses  an  algebraic  procedure  to  generate  a  finite 
difference  grid  in  the  interior,  given  the  distribution  of  mesh  points  along 
the  boundaries.  The  procedure  follows  the  transformation  credited  to  Coons 
(Reference  6)  and  described  by  Cook  (Reference  7).  The  transformation 
achieves  uniform  spacing  of  the  mesh  grid  in  the  computation  (x-y)  space,  but 
in  general  results  in  a  nonuniform,  nonorthogonal  mesh  in  the  physical  (z-r) 
space. 

There  were  a  number  of  ways  in  which  the  mesh  generator  used  in  FCT 
could  be  extended  to  include  the  chamber.  The  options  considered  for  this 
task  are  shown  in  Figure  1  along  with  the  relative  advantages  and 
disadvantages  of  each.  Generally,  extension  of  the  mesh  into  the  chamber 
results  in  poor  orientation  in  some  part  of  the  flowfield  since  the  algebraic 
generator  does  not  control  the  direction  of  the  mesh  lines,  except  those 
coincident  with  the  boundary.  An  alternative  is  to  use  an  elliptic  solver 
which  controls  both  the  spacing  and  the  direction  of  the  mesh  lines  at  and 
near  the  boundaries.  However,  adapting  an  existing  elliptic  solver  to  this 
problem  would  have  required  more  resources  than  were  available  in  this 
program.  Hence,  the  algebraic  mesh  generation  procedure  used  in  FCT  was 
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Figure  1.  Mesh  Generation  Options 


retained  for  CFC.  Moreover,  the  third  option  shown  in  Figure  1  was  selected 
because  it  yields  a  good  mesh  orientation  in  the  primary  region  of  interest, 
i.e.,  in  the  vicinity  of  the  nozzle  wall. 

The  chamber  geometry  is  input  to  CFC  by  specifying  the  radial  and 
axial  coordinates  of  various  points  along  the  boundary.  There  are 
provisions  in  CFC  to  indicate  the  presence  of  a  front  end  closure  so  that 
the  boundary  conditions  there  are  treated  differently.  The  algebraic  mesh 
generator  cannot  handle  submerged  geometries.  Hence,  a  computational 
boundary  must  be  drawn  spanning  the  submerged  region.  The  code  normally 
distributes  mesh  points  uniformly  with  the  wetted  length  along  the  chamber 
and  grain  boundary.  If  the  presence  of  a  front  end  closure  or  a  submerged 
region  inflow  boundary  is  indicated,  the  code  allows  the  specification  of  a 
minimum  number  of  mesh  points  along  the  boundaries  of  each  of  these 
regions. 

The  patch  line  separating  Regions  I  and  II  (cf..  Reference  4)  is  also 
input  by  the  user  in  CFC  by  specifying  the  coordinates  along  the  line. 

This  gives  the  user  some  control  over  the  mesh  generation  in  the  interior 
region  of  the  chamber  and  nozzle.  Except  for  the  new  features  noted  above, 
the  mesh  generation  procedure  in  the  new  code  remains  unchanged  from  that  in 
in  FCT. 

As  in  FCT,  only  axisymmetric  geometries  are  allowed  by  the  new  code. 
However,  nonaxisymmetric  configurations  such  as  star  shaped  or  axially 
slotted  grains  can  be  handled  by  approximating  them  with  equivalent 
axisymmetric  geometries.  The  procedure  to  determine  the  equivalent 
geometries  is  explained  in  detail  in  Section  5. 


2.2  BOUNDARY  CONDITIONS  IN  THE  CHAMBER 


The  gas  flow  boundary  calculations  in  CFC  have  been  modified  to  handle 
the  mass  injection  conditions  at  the  grain  and  other  inflow  surfaces.  The 
code  can  also  account  for  the  conditions  at  a  chamber  front  end  closure  with 
no  mass  injection.  These  two  types  of  boundary  conditions  are  described 
below. 

Inflow  Boundary  Conditions 

The  inflow  boundary  conditions  are  specified  in  CFC  by  inputting  a 
grain  burn-rate  law  or  mass  flux  distribution  along  the  inflow  surface. 

The  burn- rate  law  used  in  CFC  has  the  form: 

m  =  cp^  (1) 


where 


m  is  the  gas*  mass  flux  across  the  boundary 
p  is  the  static  pressure 

and 

b,  and  c  are  constants  for  a  given  location 

The  values  of  b  and  c  may  vary  with  location  along  the  boundary.  The 
burn-rate  law  can  also  be  used  to  specify  the  mass  flux  along  any  portion  of 
the  boundary  by  setting  b  =  0  in  that  region.  The  mass  flux  distribution  is 
then  input  through  c. 


*The  input  procedure  in  CFC  calls  for  the  burn-rate  law  to  be  specified  in 
terms  of  total  mass  flux  (gas  and  particles).  Since  the  particle  to  gas  mass 
ratio  is  also  input,  the  code  internally  adjusts  the  constant,  c,  to  account 
for  the  proportion  of  gas  flow  at  the  inflow  boundary. 


The  burn-rate  law  replaces  the  isentropic  relation  for  pressure  used  in 
FCT  at  the  inflow  boundary.  The  five  equations  used  to  solve  for  the  flow 
conditions  at  the  inflow  boundary  are: 

•  Compatibility  equation  along  the  negative  branch  of  the 
characteristic  curve  in  the  x-t  plane* 

•  Energy  equation  (in  the  form  of  an  adiabatic  flow  relation  for 
temperature) 

•  Equation  of  state 

•  Gas  inflow  direction 

•  Burn-rate  law 

The  form  of  the  first  four  equations  remains  unchanged  from  that  in 
FCT.  The  form  of  the  final  equation  is  derived  below.  Starting  from  the 
general  mass  flux  relation: 


m  =  pV  •  n 


we  have 


m  =  p  2ffr(udr  -  vdz) 

2irr  y4r^  +  dz^ 

where  p  is  the  gas  density,  n  is  a  unit  vector  normal  to  the  surface,  and  V 
is  the  gas  velocity  vector  with  components  u  and  v  in  the  axial  (z)  and 
radial  (r)  directions,  respectively.  A  sketch  showing  the  velocity 


*x  is  a  transformed  coordinate  which  is  zero  along  the  inflow  plane  and  one  at 
the  outflow  plane;  t  refers  to  time  (see  Reference  4). 
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components  and  defining  the  wall  angle,  0,  and  gas  flow  angle,  «,  is  shown 
be1 ow: 


Combining  Equations  (1)  and  (2)  with  the  relation  between  velocity 
components  and  the  perfect  gas  law,  we  have: 


u  =  -^  R.  (3a) 

u  ^ 


V  =  -^  R„  (3b) 

V  ^ 


where 


Ay  =  sin  8  (1  -  tan  a  cot  8)  (4a) 


Ay  =  cos  8  (cot  a  tan  8-1) 


(4b) 


Rg  is  the  gas  constant  and  T  is  the  gas  temperature. 
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The  velocity  components  given  by  Equations  {3a)  and  (3b)  satisfy  the 
relation: 

V  =  u  tan  a  (5) 

Note  that  in  CFC,  in  the  absence  of  user  input  specification  to  the  contrary, 
the  gas  flow  angle  is  assumed  perpendicular  to  the  inflow  surface. 

The  two-step  procedure  used  in  FCT  to  calculate  the  flow  conditions  at 
the  boundaries  is  maintained  in  CFC.  Within  each  step  the  five  governing 
equations  mentioned  earlier  are  solved  by  iteration  for  u,  v,  p,  T  and  p.  The 
convergence  criterion  for  the  iteration  calculations  is  based  on  pressure. 
Convergence  occurs  usually  within  two  to  three  iterations. 

No  Mass-Injection  Boundary  Condition 

A  surface  with  no  mass  injection  in  the  chamber  cannot  be  reduced  to  a 
special  case  of  the  grain  burn-rate  law  condition  described  above.  Since  such 
boundary  conditions  are  often  encountered  at  the  front  end  of  the  chamber,  CFC 
is  written  to  handle  that  region  separately,  if  required.  If  the  user  signals 
the  existence  of  a  zero  mass-injection  front  end  region  in  the  ch. mber 
geometry  input,  a  separate  set  of  equations  are  solved  by  the  code.  Five 
equations  are  again  required  to  solve  for  the  five  unknowns  (u,  v,  p,  T  and 
p),  except  at  the  center  (axis)  point  which  is  treated  as  a  stagnation  point 
with  u  =  V  =  0.  At  all  points  on  the  front  end  closure,  except  for  the 
stagnation  point,  the  five  equations  used  to  solve  for  the  flow  conditions 
are: 

1,  Compatibility  equation  along  the  negative  branch  of  the 
characteristic  curve  in  the  x-t  plane 

2.  Energy  equation  (along  a  streamline) 
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3.  Equation  of  state 

4.  Flow-surface  tangency  (slip)  condition 

5.  Compatibility  equation  along  the  negative  branch  of  the 
characteristic  curve  in  the  y-t  plane* 

The  first  and  third  conditions  remain  unchanged  from  before.  The 
fourth  condition  is  satisfied  by  making  the  flow  angle  equal  to  the  wall 
angle.  Hence, 

V  =  u  tan  B  (6) 

The  remaining  two  conditions  are  obtained  from  the  characteristic 
compatibility  equations  in  the  y-t  plane.  The  fifth  condition  is  obtained 
along  the  negative  branch  of  the  characteristic  curve,  namely: 

^  =  V  -  a5  (7) 


The  negative  (left  running)  characteristic  is  used  since  it  is  assumed  that 
the  flow  direction  is  from  the  grain  to  the  axis.  The  compatibility  equation 
along  that  characteristic  is  given  by: 


^  .  £i  Y  ^  -  £i  V  ^ 
dt  6  ^2  dt  6  ^r  dt 


£3 

5 


y,R2 


6 


yr'^3  ^  >^4 


(8) 


where 


v  =  u  y2  +  V  yp  (9) 


*y  is  the  transformed  coordinate  which  is  zero  along  the  axis  and  one  along 
the  nozzle  wall  (see  Reference  4). 
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and 


The  subscripts  2  and  r  refer  to  the  partial  derivatives  in  those  directions. 
The  term  a  is  the  local  sound  velocity.  The  variables  Rl^  R2.  R3  and  R4  in 
Equation  (8)  are  the  inhomogeneous  terms  in  the  continuity,  x-direction 
momentum,  y-direction  momentum  and  energy  equations,  respectively,  written  for 
the  y-t  reference  plane.  The  Equations  (7)  tiirough  (10)  are  from  Reference  4 
and  can  be  derived  using  the  methods  of  Reference  8. 

Equation  (8)  is  written  in  partial  differential  form*  as: 

(Pt  *  Py)  -  f-  H  ‘^y)  -  T  "y^  = 

where  x_  =  dy/dt,  given  by  Equation  (7),  and  R_  is  the  right  hand  side  of 
Equation  (8).  The  negative  subscripts  signify  that  these  terms  are  written 
for  the  negative  characteristic  branch. 

Casting  Equation  (11)  in  one-sided  difference  form,  we  have: 


Note  that  this  procedure  represents  a  departure  from  the  method  used  to 
transform  the  corresponding  equation  in  the  x-t  plane,  where  the  ordinary 
differential  of  the  compatibility  equation  along  the  characteristic  was  cast 
in  difference  form. 
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where  the  subscripts  1,  3  and  6  refer  to  the  points  shown  in  the  sketch 
below: 


The  point  3  is  at  the  new  time  where  the  flow  conditions  are  to  be 
evaluated.  Points  1  and  6  lie  in  the  same  y-t  plane  as  point  3,  but  at  an 
earlier  time  where  the  conditions  are  known.  Information  about  these 
conditions  are  transmitted  to  point  3  via  the  characteristic  in  the  y-t  plane 
shown.  The  situation  is  similar  to  that  in  the  x-t  plane  where  information  is 
transmitted  from  point  2  to  point  3  via  the  indicated  characteristic. 

However,  since  the  formulation  of  the  difference  equation  in  the  y-t  plane  was 
obtained  from  the  partial  differential  form,  the  conditions  at  point  6'  do  not 
enter  the  equation  explicitly.  Hence,  the  precise  location  of  the 
intersection  of  the  characteristic  curve  in  the  y-t  plane  with  the  y-axis, 
i.e.,  point  6',  does  not  have  to  be  evaluated. 

The  remaining  condition  (condition  2)  is  obtained  from  the  energy 
equation  along  the  streamline,  that  is: 


^  . 

dt 


,2  ^  = 


dt 


(13) 
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and 


» ■  ^  =  V  :i4) 

Again,  writing  Equation  (13)  in  partial  differential  form,  we  have: 

(Pt  ^Py)  -  ^  (pt  ""  ^4 

which  when  written  in  one-sided  difference  form  about  point  1  in  the  y-t 
plane,  gives: 


Equations  (6),  (12)  and  (16)  were  incorporated  into  the  new  code.  The 
equations  for  the  other  two  conditions  (compatibility  e^jiatior  the  x-t 
plane,  and  the  equation  of  state)  remained  unchanged  from  the  base  code.  This 
set  of  equations  are  solved  iteratively  at  the  points  along  the  boundary  where 
the  front  end  closure  is  specified.  The  convergence  criterion  is  again  based 
on  pressure. 

Finally,  at  the  center  of  the  front  end  closure,  stagnation  point 
conditions  are  imposed.  Hence  u  and  v  are  set  to  zero  and  Equations  (6)  and 
(12)  are  not  required.  The  remaining  three  equations  are  solved  iteratively 
for  pressure,  density  and  temperature,  with  the  convergence  based  on  pressure 
as  before. 
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2.3 


PARTICLE  TRAJECTORY  CALCULATIONS 


Particle  trajectory  calculations  were  modified  to  calculate  conditions 
resulting  from  trajectory  cross-overs  and  wall  impact.  Changes  were  also  made 
to  the  manner  in  which  the  particle  size  distribution  is  calculated.  These 
changes  are  discussed  below. 

Particle  Size  Distribution 

The  base  (FCT)  code  has  a  built-in  empirical  procedure  for  calculating 
the  size  distribution  of  aluminum  oxide  particles  in  the  nozzle.  However,  for 
impact  calculations  in  the  entrance  region,  the  large  aluminum  particles 
ejected  from  the  grain  are  more  important  than  the  small  aluminum  oxide 
particles.  The  large  particles  are  more  likely  to  impact  the  nozzle  walls 
since  they  cannot  follow  the  gas  stream  very  closely  around  the  nozzle 
contour.  Additionally,  larger  particles  cause  greater  mass  removal.  Hence,  a 
procedure  to  calculate  the  size  distribution  of  aluminum  particles  in  the 
chamber  was  incorporated  into  the  new  (CFC)  code. 

The  procedure  used  was  developed  by  Hermsen  (Reference  4)  for  use  in 
the  combustion  efficiency  calculations  of  the  Solid  Performance  Program  (SPP). 
The  procedure  is  not  accessible  from  the  FCT  module  of  SPP.  Hence,  it  was 
added  to  the  transonic  calculation  procedure  in  CFC.  The  average  size  of 
aluminum  particles  coming  off  the  propellant  grain  is  correlated  with  the 
grain  burnrate  and  the  percentage  of  ammonium  perchlorate  in  the  propellant. 
The  distribution  of  particle  sizes  about  the  mean  is  assumed  to  be  log-normal, 
with  an  empirically  determined  value  of  the  standard  deviation. 

Presently,  the  Chamber  Flowfield  Code  does  not  account  for  changes  in 
the  aluminum  particle  size  distribution  due  to  combustion  or  breakup.  Since 
most  of  the  breakup  occurs  in  the  throat  region  and  since  rough  calculations 
indicate  that  there  are  still  large  unburnt  particles  at  the  entrance  to  the 
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nozzle,  this  approximation  is  reasonable  for  obtaining  estimates  of  the  impact 
condition  in  the  nozzle  entrance  region. 

Trajectory  Crossovers 

The  two-phase  flow  in  FCT  is  calculated  using  a  hybrid  procedure  which 
appears  to  be  quasi-steady  but  which  actually  yields  a  fully  coupled  solution. 
The  gas  phase  equations  are  integrated  in  a  time  dependent  form,  while  the 
particle  phase  equations  are  integrated  in  a  steady  state  form.  The  gas  phase 
calculations  are  allowed  to  relax  for  a  given  number  of  time  steps  before  a 
particle  phase  calculation  is  performed.  The  gas  calculations  define  the 
flowfield  for  determining  the  aerodynamic  forces  on  the  particles.  The 
particle  calculations  in  turn  determine  the  gas-particle  coupling  terms  used 
in  the  next  set  of  gas  phase  calculation  .  The  particle  phase  calculations 
are  performed  by  tracing  the  trajectories  of  individual  particles  from 
selected  points  on  the  inflow  boundary  using  a  force-momentum  balance  on  the 
particles.  The  particle  trajectories  are  then  treated  as  particle 
streamlines.  The  spatial  density  of  particles  is  calculated  knowing  the 
particle  mass  flowrate,  the  average  velocity,  and  the  spacing  between  adjacent 
streamlines.  The  spatial  density  is  important  in  calculating  the  coupling 
terms  and  impact  conditions. 

This  hybrid  procedure  to  calculate  the  coupled  gas-particle  flow  was 
maintained  in  the  new  (CFC)  code.  However,  since  particle  trajectories  in  CFC 
start  from  the  chamber,  they  often  intersect  each  other.  Trajectory 
crossovers  may  also  occur  if  particles  cross  the  centerline,  or  if  they  impact 
the  nozzle  wall.  In  the  FCT  code  it  is  assumed  that  these  trajectory 
crossovers  do  not  occur;  therefore,  the  new  code  required  modifications  to 
account  for  the  crossovers  so  particle  spatial  densities  could  be  calculated 
correctly  with  the  total  mass  flowrate  of  particles  being  conserved. 
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The  particle  mass  flowrate  between  adjacent  streamlines  is  determined 


at  the  inflow  boundaries  by  calculating  the  gas  flowrate  between  adjacent 
trajectory  initiation  points,  and  multiplying  by  the  overall  particle  to  gas 
mass  ratio.  This  particle  mass  flowrate  associated  with  each  pair  of 
streamlines  remains  unchanged  until  crossover  or  impact  occurs.  When 
crossover  occurs,  new  pairs  of  adjacent  streamlines  are  formed  and  the  mass 
flowrate  between  them  must  be  reassigned.  An  example  of  how  this  is 
accomplished  in  CFC  is  explained  with  reference  to  the  sketch  below; 


The  sketch  shows  the  particle  trajectories  between  two  neighboring 
axial*  stations.  The  particle  densities  are  calculated  at  the  mesh  points  in 
the  code;  hence,  the  mass  flowrate  between  the  particle  streamlines  at  each 
axial  station  must  be  known.  At  each  of  those  stations,  the  particle 

*In  the  transformed  coordinates.  Hence,  axial  stations  in  this  context  refers 
to  constant  x  mesh  lines. 
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trajectory  indentifier  (index  «,),  is  renumbered  starting  from  the  bottom.  The 
code  searches  for  changes  in  the  numbering  sequence  to  indicate  the  presence 
of  particle  crossover.  If  particle  crossover  is  detected,  the  mass  flow 
between  the  streamlines  is  reassigned.  In  going  from  one  axial  station  to  the 
next  (say  from  the  line  with  index  i-1  to  the  one  with  index  i),  it  is  assumed 
that  the  particles  are  confined  between  initially  adjacent  streamlines.  This 
is  shown  schematically  below  for  each  of  the  three  pairs  of  streamlines  of  the 
earlier  sketch. 


Note  that  the  region  between  streamlines  6  and  7  at  i  gets 
contributions  from  all  the  three  flow  areas  shown  at  i-1.  Also,  the  particle 
mass  flowrate  confined  between  streamlines  5  and  6  at  i-1  is  distributed 
between  the  pairs  5-6  and  6-7  at  i .  A  similar  situation  occurs  for  the 
streamline  pair  7-8  at  i-1.  In  some  cases,  when  a  streamline  intersects  more 
than  one  other  streamline,  the  mass  between  one  pair  of  streamlines  at  one 
axial  station  may  be  spread  over  a  number  of  pairs  at  the  next  axial  station. 
The  distribution  of  mass  at  the  new  location  is  based  on  the  assumption  that 
the  particle  spatial  density  between  adjacent  streamlines  is  uniform  along  a 
line  of  constant  x.  The  mass  flowrate  between  two  streamlines.  Amp,  is 
therefore  proportional  to  the  dot  product  of  the  average  particle  velocity 
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vector,  Vp,  and  the  area  vector,  aA,  at  the  new  axial  station.  For  example, 
to  determine  the  flowrate  between  streamlines  5  and  6  at  i ,  we  have; 


%{5-6)i 


^^h5-6)i 


(17) 


al  so 


%(5-6)i-l  "  ^"’p(5-6)i  %(6-7)i 


and 


(18) 


%(6-7)i 


‘  ^^^(6-7)i 


(19) 


Hence, 


it 


AA) 


%(5-6)i  =  %(5-6)i.l-^ 


‘^^^(5-6)i 


iK  •' 


AA) 


(6-7  )i 


(20) 


A  similar  expression  can  be  derived  for  the  flowrate  between  streamlines  7  and 
8  at  i .  However,  the  flowrate  between  streamlines  6  and  7  at  i  must  account 
for  incoming  particles  from  three  flow  areas  at  i-1.  In  that  case,  we  have: 
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Note  that  this  procedure  conserves  the  particle  mass  flow  through  the  nozzle. 
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The  procedure  for  reassigning  mass  flow  between  particle  streamlines  is 
executed  within  the  code  whenever  particle  crossover  is  detected  at  or  just 
prior  to  a  given  axial  station.  When  no  crossover  occurs,  the  mass  flow 
between  adjacent  streamlines  remains  unchanged  between  subsequent  x  mesh 
locations. 

A  different  type  of  a  crossover  occurs  when  a  particle  trajectory 
intersects  the  centerline.  Since  there  is  always  a  particle  streamline  at  the 
centerline,  there  is  a  trajectory  crossover,  but  this  particular  case  is 
handled  differently  in  the  code.  When  a  particle  approaches  the  centerline, 
the  code  reflects  it  off  the  axis  in  real  space,  as  shown  in  the  sketch 
below. 


When  streamline  2  reaches  the  axis  at  C  it  is  reflected  to  point  2'  on 
the  ith  X  mesh  line.  In  the  real  flow  case  this  corresponds  to  the  streamline 
at  2  passing  down  through  the  centerline  at  C,  and  another  streamline  from  the 
opposite  side  emerging  from  the  centerline  at  the  same  location  (synmetric 
flow  assumption).  However,  in  calculating  the  particle  spatial  density  at  the 
X  mesh  lines,  the  code  assumes  the  particle  streamline  goes  directly  from  2  at 
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i-1  to  2'  at  i,  as  shown  by  the  dashed  line.  Hence,  there  is  no  need  to 
reassign  the  mass  flowrate  between  the  streamlines. 

Note  that  although  the  spatial  density  calculations  make  the  assumption 
that  the  particle  streamlines  goes  directly  from  2  to  2',  the  particle 
trajectory  calculations  do  not.  That  is,  the  portion  of  the  trajectory  2 
between  i  and  i+1  (the  straight  line  2' -2”)  is  an  extension  from  the  solid 
line  C-2',  not  the  dashed  line  2-2'. 

Wall  Impact 

The  base  (FCT)  code  does  not  allow  particles  to  impact  the  wall.  When 
particles  get  very  close  to  the  wall  they  are  forced  to  move  parallel  to  the 
wall.  Since  wall  impact  prediction  is  the  primary  objective  of  this 
analytical  procedure,  the  restriction  against  wall  impact  was  removed  in  the 
new  (CFC)  code.  Further,  the  mass  deposited  on  the  walls  and  the  conditions 
of  impact  are  calculated. 

The  mass  of  particles  reaching  the  wall  are  calculated  using  a 
procedure  similar  to  that  used  in  distributing  the  mass  flow  after  trajectory 
crossover.  The  location  of  impact  is  calculated  and  stored  by  the  code.  When 
the  code  begins  to  calculate  mass  flow  between  streamline  and  particle  spatial 
densities,  it  checks  for  wall  impact  (in  addition  to  trajectory  crossover) 
from  one  axial  station  to  another.  If  impact  is  detected,  the  particle  mass 
flowrate  reaching  the  wall  is  calculated.  It  is  assumed  that  the  mass  of 
particles  reaching  the  wall  is  removed  from  the  flow. 

The  procedure  for  calculating  the  mass  of  particles  reaching  the  wall 
is  described  with  reference  to  the  example  shown  in  the  sketch  below.  The 
particle  streamline  10  impacts  the  nozzle  wall  at  point  P  between  the  axial 
stations  i-1  and  i.  The  mass  flowrate  between  the  streamlines  9  and  10  at  i-1 
is  divided  into  two  parts:  some  of  the  mass  is  deposited  on  the  wall  between 
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P  and  10'  and  the  remainder  passes  through  the  points  9'  and  10'.  (The 
point  10'  is  where  the  axial  station  i  intersects  the  nozzle  wall).  The 
distribution  of  mass  is  again  taken  proportional  to  the  values  of  the  vector 
dot  product  Vp  .  dA  between  P  and  10'  and  between  9'  and  10'.  The  vector  dot 
product  is  based  on  the  average  velocity  in  the  region  of  interest.  The 
average  velocity  is  obtained  by  taking  the  arithmetic  mean  of  the  velocity 
components  at  the  regions'  extremities.  In  this  case,  the  velocity  at 
point  10'  is  not  known.  It  is  assumed  to  be  the  same  as  point  P  since  the  two 
points  are  generally  located  very  close  to  one  another. 

The  procedure  is  then  repeated  at  the  next  axial  station.  The  mass 
flow  through  9'  and  10'  at  i  is  partly  deposited  at  the  wall  between  the 
points  10'  and  10",  while  the  rest  flows  through  the  region  between  the 
points  9"  and  10".  Note  that  at  each  of  the  downstream  x-mesh  wall 
points  10',  10",  etc.,  there  is  a  remnant  of  the  original  streamline  10.  The 
procedure  continues  at  subsequent  axial  stations  until  another  streamline 
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impacts  the  wall  or  the  remnant  streamline  breaks  free  from  the  wall.  The 
latter  occurs  when  the  wall  curves  away  sufficiently  that  the  remnant  velocity 
vector  is  no  longer  pointing  into  the  nozzle  wall.  Once  a  streamline  breaks 
free  from  the  wall,  it  is  treated  as  any  other  streamline  in  the  flow,  with 
the  mass  between  it  and  the  adjacent  streamline  reduced  to  account  for  mass 
deposition  upstream  of  the  breakaway  point. 

As  mentioned  earlier,  the  code  calculates  impact  conditions  as  well  as 
the  impact  location.  The  quantities  required  for  calculating  erosion  due  to 
particulate  impact,  such  as  particle  velocity,  mass  flux  and  impact  angle,  are 
calculated  and  printed  in  the  output  for  each  particle  size  group.  The 
particle  size  group  gives  the  remaining  parameter  required  for  calculating 
erosion,  that  is,  particle  size.  The  impact  conditions  are  then  used  as  input 
to  the  ablation-erosion  code  which  is  discussed  in  Part  2. 


p.  3  (»  B 
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SECTION  3 


CFC  OVERALL  STRUCTURE 


The  Chamber  Flowfield  Code  (CFC)  is  structured  as  a  stand-alone  version 
of  the  Fully  Coupled  Transonic  (FCT)  module  of  the  Solid  Performance  Program 
(SPP).  Figure  2  shows  a  comparison  between  the  overall  structures  of  SPP  and 
CFC.  The  control  module  in  CFC  serves  only  to  read  in  the  variables  in  the 
SPP  namelist  required  for  CFC,  and  call  the  TD2P  module  which  in  turn  executes 
the  modified  FCT. 

The  input  data  requirements  for  CFC  are  given  in  Table  1.  A  brief 
description  of  each  input  is  also  given,  and  sections  containing  greater 
detail  are  identified.  The  input  data  set  for  CFC  contains  fewer  items  than 
the  complete  set  for  SPP.  However,  note  that  there  is  a  new  namelist,  called 
TD3,  where  the  new  input  variables  required  for  the  modified  version  of  FCT 
are  specified. 

Since  CFC  does  not  contain  thermochemical  expansion  or  ballistics 
modules,  the  variables  transmitted  in  SPP  from  these  modules  to,  and  used  by, 
FCT  must  now  be  specified  in  the  TD2  namelist.  The  only  variable  transmitted 
from  the  ballistics  module  is  the  chamber  pressure.  This  can  be  obtained  from 
a  test  firing  pressure  trace,  if  available,  or  by  running  a  ballistics  code, 
such  as  the  BAL  module  of  SPP.  The  thermochemical  variables  can  be  obtained 
by  running  an  equilibrium  expansion  code,  such  as  the  One-Dimensional 
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Figure  2.  Comparison  of  Overall  Structures  of  SPP  and  CFC 


Table  1.  Input  Data  Set  Description 


Card  Input 

Description 

Section  No. 

T  I  T  L  E  [ 

One  or  more  title  cards 

4.0 

SPP 

SPP  directive  card 

4.0 

$  S  P  P  1 

SPP  namelist  for  nozzle  geometry 

4.1 

:  I 

and  computer  plotting 

$  E  N  D  ) 

$  T  D  2  1 

: 

T02  namelist  for  old  FCT  variables 

4.2 

send) 

$  T  D  3  ) 

TD3  namelist  for  new  FCT  variables 

5.0 

send) 

Equilibrium  (ODE)  module  of  SPP  or  the  Aerotherm  Chemical  Equilibrium  (ACE) 
code  (Reference  9).  An  example  of  how  an  ACE  expansion  is  utilized  for  CFC 
input  is  shown  in  the  discussion  of  the  sample  case  (Section  7). 
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SECTION  4 


SPP  INPUT  CHANGES  FOR  CFC 

This  section  describes  changes  to  the  existing  SPP  input  variables 
required  to  execute  CFC.  Essentially  the  input  is  set  up  as  if  FCT  alone  were 
to  be  run  on  SPP.  Of  all  the  data  required  by  the  SPP  code,  only  the  SPP 
directive  card  and  the  SPP  and  TD2  namelists  are,  therefore,  required.  Title 
cards  may  be  input,  if  desired.  The  formats  for  the  title  cards,  the  SPP 
d’  "ective  card  and  the  two  namelists  are  the  same  as  in  SPP.  However,  the 
namelist  data  requirements  are  simplified  because  only  the  modified  FCT  module 
is  executed,  and  since  specific  impulse  calculations  are  bypassed. 

4.1  SPP  NAMELIST 

The  function  of  the  SPP  namelist  in  CFC  is  to  input  the  nozzle  wall 
geometry  and  computer  plotting  variables.  Changes  to  the  namelist  are 
describ<:>J  below  under  the  same  headings  used  in  the  SPP  user's  manual 
(Reference  4). 

Modules  to  be  Executed,  Etc. 

TD2P  =  1. 

This  is  presently  the  only  acceptable  specification.  Calls  to  other 
modules  have  been  eliminated. 

ODK  Start  Option 

Not  required. 
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Missile  Trajectory 


Not  required. 

Nozzle  Geometry 

Specify  the  following: 

RSI 

Nozzle  Inlet  Geometry  (INLET  and  corresponding  geometry  variables) 
Nozzle  Exhaust  Geometry  (IWALL  and  corresponding  geometry  variables) 
Note:  Although  CFG  calculations  extend  only  up  to  the  nozzle  attachment 
point,  a  cone  or  contour  geometry  beyond  that  location  must  be  specified. 
The  geometry  beyond  the  attachment  point  may  be  fictitious  since  it  will 
not  affect  the  calculations. 

Nozzle  Throat-Erosion  Data 
Not  required. 

Combustion  Efficiency  Calculations 

Not  required.  Calculation  of  aluminum  agglomerate  sizes  using  the 
internal  empirical  correlation  is  available  as  an  option  in  the  TD3 
namel ist. 

Computer  Plotting 

FCT  plots  only.  Specify  PLTT,  PRNTPL  and  FCT  plot  variables. 

Note:  No  changes  have  been  made  to  the  SPP  plot  routines.  The  code,  as 
received,  does  not  save  plots.  Although  the  plots  are  not  saved,  both  the 
NPLTS  and  NSPLTS  vectors  must  be  specified  for  generating  plots  during  a 
run.  Set  NSPLTS  =  NPLTS. 

4.2  TD2  NAMELIST 

Changes  to  the  TD2  namelist  are  described  under  the  same  headings  as 
used  in  the  SPP  user's  manual. 
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Items  Communicated  for  ODE  and  BAL  Modules 


All  items  must  be  specified.  In  addition,  CPS,  the  solid  particle 
heat  capacity  in  ft^/sec^  °R  should  be  included  in  the  namelist. 

The  default  value  of  XMLW,  the  molecular  weight  of  the  condensed 
phase,  has  been  changed  to  that  for  aluminum  (26.98154  Ib/lbmole). 

Perfect  Gas  Option 

If  the  perfect  gas  option  is  desired  (P6FLG  is  set  equal  to  1.), 
GAMMA  and  CPG  must  be  specified. 

Particle  Data 

The  five  options  for  calculating  aluminum  oxide  particle  group 
sizes  in  SPP  are  still  in  effect  in  CFC.  However,  the  default  values  of 
some  of  the  variables  have  been  changed  to  facilitate  aluminum  agglomerate 
particle  group  size  calculations.  Options  4  and  5  in  SPP  can  also  be  used 
in  CFC  to  directly  input  aluminum  particle  sizes  and  distributions.  In 
addition,  a  sixth  option  has  been  added  to  calculate  the  mean  aluminum 
particle  size  via  an  empirical  correlation.  These  options  are  discussed 
bel ow. 

Option  4:  Input  of  DPBAR 

In  this  option  a  mean  particle  diameter,  DPBAR,  is  input.  If  the 
number  of  particle  groups  re<juested  is  more  than  one,  the  particle  group 
sizes  are  automat ica' 1 y  calculated  using  a  log  normal  distribution.  The 
number  of  particle  groups  desired  is  input  through  the  variable  NPG,  and 
the  geometric  standard  deviation  of  the  log  normal  distribution  is  input 
through  SIG.  As  in  SPP,  the  particle  mass  fractions  are  distributed 
uniformly  over  all  size  groups. 
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Option  5:  Input  of  R  and  WPWT 

In  this  option  the  size  of  the  particle  groups  are  input  through  the  R 
array,  and  the  corresponding  weight  fractions  of  each  group  are  input  through 
the  WPWT  array. 

Option  6:  Input  of  ROOT  and  ALFAAP  (in  TD3  namelist) 

This  is  the  recommended  option.  Here  the  mean  particle  diameter, 

DPBAR,  is  not  input  but  calculated  using  Hermsen's  correlation  for  mean 
aluminum  agglomerate  size.  This  correlation,  which  is  also  used  for 
combustion  efficiency  calculations  in  SPP,  relates  DPBAR  to  the  grain  burn 
rate  (ROOT)  and  to  the  mass  fraction  of  ammonium  perchlorate  in  the  propellant 
(ALFAAP)  as  follows: 

DPBAR  =  35.0/RD0T/ALFAAP 

The  particle  size  groups  are  then  calculated  using  the  log  normal 
distribution.  Hermsen  characterizes  the  aluminum  agglomerate  size 
distribution  by  a  standard  deviation  of  logjo^^  equal  to  0.2.  The  default 
value  of  SIG  was  selected  to  be  consistent  with  this  empirically  determined 
val ue.* 

Options  4  Through  6 

In  addition  to  the  variables  mentioned  above,  the  particle  material 
density  is  required.  It  can  be  in  terms  of  the  variable  SMP  if  the  material 
densities  of  all  particle  groups  can  be  considered  equal  and  invariant  with 
temperature.  Later,  an  array  variable  ROPMK  is  described  which  allows  for  the 
specification  of  density  as  a  function  of  particle  group  and  temperature. 


*SIG  is  related  to  og,  the  standard  deviation  of  logjoD.  as  follows: 
SIG  =  lo'^a 
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The  code  will  assume  a  default  value  of  147.6  Ibm/ft^  (density  of  liquid 
aluminum),  if  SMP  or  ROPMK  are  not  specified  in  the  TD2  namelist. 

The  input  variables  for  options  4  through  6  are  summarized  below; 


Option 

Item 

Description 

Units 

Default 

Value 

4,6 

NPG 

Number  of  particle  groups  to  be 
considered.  NPG  must  be  less  than 
or  equal  to  5. 

None 

3 

4  to  6 

SMP 

Particle  material  density 

1 bm/ft^ 

147.6 

4 

OPBAR 

Mean  particle  diameter 

m 

0.0 

4,6 

SIG 

Geometric  standard  deviation  for 
log  normal  particle  size  distribution 

None 

1.585 

5 

R(l) 

Array  of  particle  radii.  Maximum 
of  5.  Must  be  input  in  ascending 
orcleT^.  Add  a  zero  to  the  list  after 
the  final  (maximum)  value  to  indicate 
end  of  set. 

ft 

0.0 

5 

WPWT(l) 

Particle  weight  flow  fractions 
corresponding  to  each  of  the  above 
particle  radii.  Entries  must  add 
to  one. 

None 

6 

ROOT 

Propellant  burning  rate  (input  in 

TD3  namelist) 

in ./sec 

0.0 

6 

ALFAAP 

Mass  fraction  of  ammonium  perchlorate 
in  propellant  (input  in  TD3  namelist) 

None 

0.0 

NOTE 


If  OPBAR  is  specified  (other  than  zero),  Option  4  will  be 
executed  and  inputs  of  R(l)  or  ROOT  and  ALFAAP,  if  any,  will 
be  overridden.  Also,  an  input  of  R(l)  will  execute  Option  5 
only  if  OPBAR  is  zero,  regardless  of  ROOT  and  ALFAAP  values. 
Finally,  to  execute  Option  6,  both  ROOT  and  ALFAAP  must  be 
different  from  zero,  and  OPBAR  and  R(l)  must  both  be  zero. 
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Transonic  Analysis  Input  Data 

The  input  data  used  in  CFC  are  given  below: 


Item 

Description 

Units 

Default 

Value 

FCT 

.TRUE.  The  user  has  no  discretion  over 
this  input.  The  Advanced  Transonic 
Analysis  (ATA)  option  is  not  available 
in  CFC 

logical 

.TRUE. 

RDPMK(1,K) 

147.6 

Particle  material  density  table  for  the 
particle  group 

1 bm/ft^ 

TROMK  (1,K) 
0.0 

Temperature  table  corresponding  to  the 
particle  density  table 

°R 

NROM  (K) 

Number  of  entries  in  the  ROPMK  (1,K) 

TROMK  (1,K)  tables.  NROM  (K)  _<  10 

none 

0.0 

NSMAX 

Maximum  number  of  time  steps  to  be  taken 

none 

1,000 

NGIPI 

Number  of  gas  phase  calculation  steps 
per  set  of  particle  phase  calculations 

none 

25 

Print  Control 

Same  as  in  SPP. 

Time  Step  Control 

Same  as  in  SPP. 

Mesh  Generation  Data 

The  inlet  boundary  and  the  patch  line  demarcating  Regions  I  and  II 
are  input  in  the  T03  namelist.  The  only  variables  that  need  to  be  input 
in  the  TD2  namelist  are:  NX,  NY,  REG2,  and  ZMAXA.  The  description,  units 
and  default  values  of  these  variables  are  the  same  as  in  SPP.  It  is 
recommended  that  the  variable  REG2  be  set  to  .TRUE,  and  the  patch-line 
input  as  described  in  Section  5. 
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Inflow  Boundary  Condition  Data 

The  gas-flow  inlet  boundary  conditions  are  specified  in  terms  of  a 
grain-burn  law  or  mass-flux  distribution  in  the  TD3  namelist.  Hence,  the 
option  to  specify  any  desired  total  pressure  and  temperature  variation 
along  the  inflow  boundary  (via  the  vectors  POIN  and  TOIN)  is  no  longer 
applicable.  The  following  options  remain  operational  in  the  TD2  namelist 
for  specifying  the  gas-phase  boundary  conditions: 

ITHIN  =  flag  controlling  the  distribution  of  the  gas-phase  flow  angle 
along  the  inflow  boundaries  in  the  chamber 
=  -1  The  flow  angle  is  assumed  normal  to  the  inflow  boundary 

(default  option) 

=  0  No  longer  applicable 

=  1  The  flow  angles  are  input  via  the  THIN  array 

If  ITHIN  is  set  equal  to  1,  the  table  THIN  is  specified  in  the  same 
manner  as  in  SPP.  In  this  case,  the  vectors  POIN(J)  and  TOIN(J)  must  be 
specified  and  set  equal  to  1.0.  Note  that  in  the  default  case 
(ITHIN  =  -1)  the  gas  flow  angle  is  assumed  normal  to  the  inflow  boundary, 
which  is  a  departure  from  the  procedure  used  in  SPP.  The  vectors  POIN  and 
TOIN  need  not  be  specified  for  the  default  option. 

The  data  set  for  the  particle  phase  boundary  conditions  remains  the 
same  as  in  SPP.  The  particle  velocity  is  specified  as  a  fraction  of  the 
gas  phase  velocity  via  the  array  PKK.  The  default  value  for  PKK  is  1.  It 
is  recommended  that  PKK  be  set  equal  to  0.2  to  account  for  the  velocity 
lag  of  the  aluminum  agglomerate  at  the  propellant  surface. 

Remaining  TD2  Data  Sets 
Not  required. 
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SECTION  5 


NEW  INPUT  VARIABLES  FOR  CFC 


New  input  variables  required  for  the  Chamber  Flowfield  Code  are 
read  in  through  the  namelist  TD3.  The  new  variables  are  related  to 
changes  in  geometry,  inflow  boundary  conditions  and  particle  trajectory 
calculations.  The  data  sets  for  these  and  other  variables  are  described 
bel ow. 

Geometry 

The  chamber  geometry  is  input  in  CFC  by  specifying  the  boundary 
coordinates  RCHM8  and  ZCHMB,  for  the  r  and  z  directions,  respectively.  In 
general,  the  input  coordinates  describe  the  actual  chamber  and  grain 
surface.  However,  in  some  cases  they  may  define  hypothetical 
computational  inflow  boundaries.  For  example,  an  axisymmetric 
approximation  to  the  real  boundaries  must  be  employed  for  star-shaped, 
slotted,  or  other  nonaxi symmetric  grain  configurations.  The  recommended 
procedure  is  as  follows:  at  a  given  axial  location,  ZCHMB,  an  equivalent 
axisymmetric  (circular)  grain  surface  is  defined  with  a  port  area  equal  to 
the  actual  axial  flow  area.  The  radius  of  the  equivalent  circular 
boundary  is  the  radial  coordinate,  RCHMB,  corresponding  to  ZCHMB. 

Figure  3  explains  the  details  of  the  procedure.  The  mass  injection  rate 
must  also  be  adjusted  to  account  for  the  difference  in  surface  areas 
between  the  real  and  computational  boundaries.  A  perimeter  adjustment 
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factor,  FPERCH,  is  used  for  this  purpose.  As  explained  in  Figure  3,  FPERCH  is 
the  ratio  of  the  real  surface  area  to  the  computational  boundary  area  at  a 
given  axial  location. 

The  factor  FPERCH  may  also  be  used  in  other  regions  of  the  chamber 
where  the  computational  boundary  does  not  correspond  to  the  real  boundary. 

For  example,  the  current  mesh  generator  has  difficulty  in  mapping  the 
submerged  region.  A  computational  boundary  must  therefore  be  drawn  connecting 
the  nosecap  tangent  point  to  some  point  on  the  grain  surface,  as  shown  by  line 
DE  in  Figure  4.  It  is  recommended  that  the  line  DE  should  not  be  at  an  angle 
greater  than  60°  to  the  nozzle  axis.  The  factor  FPERCH  is  then  used  to  adjust 
the  mass  flow  across  the  computational  boundary  DE  so  that  it  approximates  the 
mass  injected  by  the  grain  in  the  submerged  region  DF. 

The  chamber  coordinates  RCHMB  and  ZCHMB  are  input  via  two 
one-dimensional  tables,  similar  to  the  nozzle  wall  coordinates  specification. 
Likewise,  the  coordinates  are  normalized  by  throat  radius,  and  the  origin  is 
at  the  center  of  the  throat.  The  coordinates  are  specified  starting  from  the 
front-end  center  of  th*.  chamber  (point  A  in  Figure  4).  The  last  entry  in  the 
tables  corresponds  to  the  point  where  the  chamber  geometry  meets  the  nozzle 
wall  (point  E  in  Figure  4).  The  total  number  of  points  input  is  specified  by 
the  variable  ICHMB. 

In  many  motor  designs,  the  front  end  of  the  chamber  has  an  impermeable 
boundary  such  as  a  solid  wall;  e.g.,  line  AB  in  Figure  4.  In  that  case  the 
extent  of  the  no- injection  region  must  be  specified  so  that  the  code  treats 
that  portion  of  the  boundary  separately,  using  slip-flow*  instead  of 


*The  point  at  the  center  of  the  front-end  closure  is  treated  as  a 
stagnation  point.  Therefore,  in  the  chamber  geometry  specification, 
the  contour  next  to  point  A  in  Figure  4  should  be  normal  to  the  wall. 
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Figure  4.  Chamber  and  Grain  Configuration  Schematic 


mass-injection  boundary  conditions.  This  can  be  done  by  including  the 
coordinates  of  the  point  demarcating  the  no- injection  and  mass-injection 
boundaries  in  the  RCHMB  and  ZCHMB  tables  (point  B  in  Figure  4).  The  point 
is  then  identified  by  specifying  its  position  on  the  table  via  the  variable 
ICHMBl.  If  a  value  of  ICHMBl  is  specified  (other  than  1),  the  code  will 
automatically  locate  a  mesh  point  at  the  demarcation  point.  The  minimum 
number  of  mesh  points  in  the  no-injection  region  can  also  be  input,  if 
desired,  by  specifying  the  variable  NYMINl. 

The  extent  of  the  imaginary  boundary  spanning  the  submerged  region 
(line  DE  in  Figure  4)  can  also  be  specified  in  a  similar  manner.  The 
coordinates  of  point  D  are  identified  by  specifying  their  array  subscript, 
ICHMB2.  A  mesh  point  will  then  be  automatically  placed  at  D.  The  minimum 
number  of  mesh  points  along  OE  can  also  be  input  by  specifying  NYMIN2. 

If  Region  II  mesh  generation  is  desired  (REG2  =  .TRUE,  in  the  TD2 
namelist),  the  user  must  also  input  normalized  coordinates  of  the  patch  line 
separating  Regions  I  and  II.  The  patch  coordinates  are  input  via  two 
one-dimentional  arrays:  RPATCH  for  the  r  coordinate,  and  ZPATCH  for  the  z 
coordinate.  The  total  number  of  coordinate  pairs  input  is  specified  by  the 
variable  IPATCH. 

Line  XY  in  Figure  4  shows  a  typical  patch  line.  The  location  of  X  and 
Y  is  such  that  the  ratios  of  the  lengths,  OX/OA  and  PY/PE,  are  approximately 
equal.  This  ensures  that  the  mesh  spacing  along  the  axis  and  the  nozzle  wall 
is  approximately  uniform  in  Regions  I  and  II.  Generally  X  and  Y  should  be 
located  such  that  the  length  ratios,  OX/OA  and  PY/PE,  are  between  one-quarter 
and  one-half.  Once  the  points  X  and  Y  are  located,  a  smooth  curve  is  drawn 
joining  the  two  points.  The  curve  should  be  roughly  normal  to  the  axis  and 
the  wall  at  points  X  and  Y,  respectively. 
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The 

table  below  summarizes  the  geometry  variables  input  in  the  TD3 

name! i st . 

Item 

Description 

Units  Default  Value 

RCHMB 

table  of  normalized  chamber  (front- 
closure-grain-infl ow  boundaries) 
radial  coordiates 

non- 

dimen¬ 

sional 

ZCHMB 

table  of  normalized  chamber  axial 
coordinates 

non- 

dimen¬ 

sional 

ICHMB 

number  of  entries  in  the  RCHMB, 

ZCHMB  tables.  ICHMB  _<50 

none 

ICHMBl 

subscript  of  the  chamber  coordinates 
at  the  junction  of  a  no-injection 
front  end  closure  and  a  mass 
injection  boundary 

none 

% 

ICHMB2 

subscript  of  the  chamber  coordinates 
at  the  start  of  the  imaginary  inflow 
boundary  across  the  submerged  region 

none  ICHMB 

NYMINl 

minimum  number  of  mesh  points  along 
the  front-end  closure  boundary 
(NYMINl  >2,  if  ICHMBl  >1) 

none  3 

NYMIN2 

minimum  number  of  mesh  points  along 
the  imaginary  boundary  across  the 
submerged  region  (NYMIN2  >2,  if 

ICHMB2  <  ICHMB) 

none  2 

RPATCH 

table  of  normalized  patch  (line 
separating  Regions  I  and  11)  radial 
coordinates 

non- 

dimen¬ 

sional 

ZPATCH 

table  of  normalized  patch  axial 
coordinates 

non- 

dimen¬ 

sional 

I  PATCH 

number  of  entries  in  the  RPATCH, 
ZPATCH  tables.  I PATCH  _<50 

none 

t 
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Inflow  Boundary  Conditions 

The  inflow  boundary  conditions  are  specified  in  the  form  of  a  mass 
burn-rate  law  or  mass  flux  distribution  along  a  grain  or  other  inflow  surface. 
The  burn-rate  law  has  the  expression: 

m  =  cp^ 

where  m  is  the  mass  flux  in  Ibm/ft^-sec 
p  is  the  pressure  in  psia 
and  b,  c  are  constants  at  any  given  location. 

The  mass  flux  specification  option  is  treated  as  a  special  case  of  the 
burn-rate  law  by  setting  the  constant  b  equal  to  zero. 

Since  the  constants  b  and  c  may  vary  along  the  grain  and  other  inflow 
boundaries,  they  are  input  in  tabular  form  via  two  one-dimensional  arrays: 
8EXPCH  for  b  and  CBRNCH  for  c.  The  values  of  the  constants  are  input  in 
sequence  corresponding  to  each  point  in  the  chamber  coordinate  tables;  that 
is,  BEXPCH  (I)  and  CBRNCH  (Z)  are  the  values  of  b  and  c  at  the  point  RCHMB 
(I),  ZCHMB  (I).  The  values  of  the  constants  between  two  input  points  are 
calculated  within  the  code  by  linear  interpolation  along  the  length  of  the 
chamber  boundary. 

The  factor  FPERCH  (see  the  discussion  under  geometry)  is  also  input  in 
a  one-dimensional  array.  The  sequence  of  values  in  the  FPERCH  table 
corresponds  to  the  sequence  of  points  in  the  chamber  coordinate  tables.  The 
code  adjusts  the  mass  injection  rate,  as  calculated  by  the  burn-rate  law,  by 
multiplying  it  by  the  corresponding  value  of  FPERCH. 
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The  inflow  boundary  conditions  input  in  the  TD3  namelist  are  summarized 


bel ow. 

Item 

Description 

Units 

Defaul  t 

Val 

BEXPCH 

exponent  b  used  in  the  grain 
burn-rate  law 

none 

50  * 

0.0 

CBRNCH 

constant  c  used  in  the  grain 
burn-rate  law 

bm/ft^  -sec 
(psia)'’ 

50  * 

1.0 

FPERCH 

ratio  of  actual  to  input  surface 
areas 

none 

50  * 

1.0 

NOTE 

The  above  three  quantities  are  input  in  tabular  form,  with  the 
values  corresponding  sequentially  to  the  points  in  the  chamber 
geometry  coordinate  tables. 

Particle  Size  Calculation 

If  option  6  is  selected  for  calculating  the  particle  size  groups,  the 
values  of  ROOT  and  ALFAAP  should  be  input  (different  from  zero)  in  the  TD3 
namelist.  See  Section  3  for  details. 

Particle  Streamline  Starting  Location  Distribution 

The  distribution  of  the  starting  points  of  particle  streamlines  along 
the  inflow  boundary  can  be  controlled  by  the  user.  The  particle  streamline 
starting  locations  are  based  on  their  particle  mass  flow  fractions.  Each 
streamline,  L,  has  a  given  initial  mass  flow  fraction,  WPSCOL(L).  The  Lth 
streamline  starts  from  a  location  such  that  the  fraction  of  th‘;  ‘'.al  particle 
mass  flow  injected  between  it  and  the  axis,  along  the  inflow  boundary,  is 
equal  to  SPSCDL(L).  The  first  particle  streamline  transverses  the  centerline; 
hence  WPSCOL(l)  =  0.0.  Also  the  final  streamline  is  a  limiting  streamline, 
hence  WPSCDL{NIPPG)  =  1.0,  where  NIPPG  is  the  total  number  of  streamlines  (set 
equal  to  10  in  the  code).  The  other  values  of  WPSCDL(L)  range  between  0  and 
1,  arranged  in  ascending  order  of  N. 
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The  vector  WPSCOL  may  be  directly  read  in,  if  desired,  in  the  TD3 


namelist.  Another  available  option  is  to  let  the  code  generate  the 
distribution  internally  according  to  the  relation: 

WPSCDL{L)  .  {L-1)WPEXP 

where  the  exponent  WPEXP  is  an  input  constant.  If  the  WPSCDL  array  is  read  in 
directly,  WPEXP  must  be  input  equal  to  zero. 

If  the  automatic  WPSCOL  generation  option  is  used,  the  minimum  number 
of  streamlines  starting  from  the  submerged  region  (i.e.,  the  imaginary  inflow 
boundary  spanning  that  region)  can  be  specified.  To  use  this  feature,  WPEXP 
must  be  specified  other  than  zero,  and  ICHMB2  must  be  input  to  define  the 
extent  of  the  submerged  inflow  region.  The  minimum  number  of  streamlines  is 
then  specified  via  the  variable  NPSSUB.  NPSSUB  is  always  greater  than  or 
equal  to  2. 

The  variables  used  to  define  the  particle  streamline  distribution  at 
the  inflow  boundary  are  summarized  below. 


Item 

Description 

Units 

Default  Value 

WPEXP 

exponent  used  to  automatically 
generate  the  WPSCDL  array 

none 

2.0* 

NPSSUB 

minimum  number  of  particle  stream¬ 
lines  starting  from  the  submerged 
inflow  boundary 

none 

2 

WPSCDL 

table  of  particle  mass  flow 
fractions  for  each  streamline  at 
the  inflow  boundary 

none 

10  *  0.0 

*When  WPEXP  is  set  equal  to  2.0,  the  code  generates  values  of  WPSCDL 
approximately  equal  to  those  used  in  the  original  FCT  module. 
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Output  Control 

If  particle  impact  occurs,  the  impact  conditions  are  automatically 
printed  out  with  the  regular  output  (at  the  output  steps  specified  by  the  NPRT 
vector  in  the  TD2  namelist).  The  user  can  also  ask  for  a  printout  detailing 
the  particle  streamline  positions  and  crossover  conditions  at  each  output  step 
by  betting  the  logical  variable  PTPART  to  .TRUE. 

The  code  also  has  provisions  to  print  a  detailed  debug  output  which 
lists  each  iteration  of  the  trajectory  integration  calculations.  The  debug 
output  is  obtained  by  setting  the  logical  variable  PTDTRJ  to  .TRUE.  The 
reader  is  cautioned  not  to  use  the  latter  option  indiscriminately  since  a 
massive  output  results. 

A  summary  of  the  output  control  variables  is  given  below. 


Item 

Description 

Units 

Default  Value 

PTPART 

=  .TRUE.,  will  result  in  a  detailed 
printout  of  particle  streamline 
positions  and  crossover  conditions 
at  each  output  step 

logical 

.FALSE. 

PTDTRJ 

=  .TRUE.,  will  result  in  a  debug 
printout  of  all  trajectory 
integration  iterations  at  each 
output  step 

1 ogical 

.FALSE. 

Miscellaneous  Variables 

There  are  a  number  of  other  variables  in  the  TD3  namelist  used  to 
facilitate  calculations.  Numerical  instabilities  in  the  gas  flowfield  are 
controlled  by  the  variable  DAMPN.  If  DAMPN  is  set  to  .TRUE.,  the  gas 
pressure,  density,  and  velocity  components  are  smoothed  using  a  fourth-order 
scheme  described  by  Kutter  (Reference  10).  The  smoothing  is  controlled  by  the 
coefficients  SMLDI  and  SMLDO  in  the  x  and  y  directions  of  the  transformed 
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coordinate  plane.  The  coefficients  range  in  value  from  0.0  (no  smoothing)  to 
0.5  (maximum  smoothing).  Values  between  0.05  and  0.3  are  recommended. 

Instabilities  in  the  gas  flowfield  are  especially  critical  during  the 
early  phase  of  the  calculation.  Stability  is  enhanced  by  allowing  the  gas 
flowfield  to  relax  for  a  sufficient  number  of  time  steps  before  particle-gas 
coupling  terms  are  introduced  in  the  calculations.  The  variable  NGIPIl  delays 
the  particle  trajectory  calculations  until  NGIPIl  gas  flowfield  time  steps 
have  been  executed.  Values  between  100  and  200  for  NGIPIl  are  recommended. 

Savings  in  computational  time  can  be  effected  by  reducing  the  number  of 
particle  trajectory  calculations.  The  number  of  gas  phase  calculation  steps 
per  set  of  particle  phase  calculations,  defined  by  the  variable  NGIPl  in  the 
T02  namelist,  has  been  increased  to  25  in  CFC,  thus  decreasing  the  total 
number  of  particle  trajectory  calculations.  These  calculations  can  be  further 
reduced  by  specifying  only  one  particle  size  group.  However,  for  impact 
calculations,  trajectories  of  more  than  one  size  group  are  of  interest.  By 
setting  the  variable  SPCUPL  to  .TRUE.,  the  time-dependent  calculations  are 
performed  with  only  one  particle  group,  while  the  calculations  at  the  final 
time  step  are  performed  with  NPG  particle  groups.  (The  variable  NPG  is 
specified  in  the  TD2  namelist.)  The  user  can  control  the  particle  group  used 
for  the  single-sized  particle  calculations  through  the  variable  KSPAR.  It  is 
recommended  that  KSPAR  be  set  to  the  particle  group  corresponding  to  the  mean 
particle  size. 

The  variable  WTROX  is  used  to  adjust  the  mass  flux  of  aluminum 
particles  impacting  the  wall.  The  code  currently  does  not  account  for  mass 
transfer  beween  the  phases.  Since  it  is  important  to  obtain  the  correct 
choking  conditions ,  and  since  the  particle  phase  at  the  throat  is  comprised 
mostly  of  aluminum  oxide,  the  variable  WPWGT*  should  be  based  on  the  aluminum 
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oxide-to-gas  mass  ratio  at  the  throat.  However,  since  impact  by  aluminum 
agglomerate  is  of  interest  in  the  entrance  region,  the  mass  of  particles 
impacting  the  wall  must  be  adjusted  by  the  AI/AI2O3  mass  ratio.  This  ratio  is 
specified  in  the  TD3  namelist  by  the  variable  WTROX. 

A  summary  of  variables  is  given  below. 


Item 

Description 

Units 

Default  Value 

DAMPN 

=  .TRUE.,  will  result  in  numerical 

smoothing  of  gas  flowfield  variables 

logical 

.FALSE. 

SMLOI 

smoothing  coefficient  in  x  direction 
of  transformed  coordinate  plane 

none 

0.2 

SMLDJ 

smoothing  coefficient  in  y  direction 
of  transformed  coordinate  plane 

none 

0.2 

NGIPIl 

number  of  gas  phase  calculation 
steps  before  first  particle 
phase  calculation 

none 

0 

SPCUDL 

=  .TRUE.,  will  result  in  one  particle 
size  group  being  used  in  all 
calculations,  except  for  the 
final  step 

none 

.FALSE. 

KSPAR 

the  particle  size  group  to  be  used 
in  above  calculation 

1  _<  USPAR  _<  NPG 

none 

1 

WTROX 

metal -to-metal  oxide  mass  ratio 

WTROX  =  0.5292  for  aluminum 

none 

1.0 

♦WPWGT  is  the  particle-to-gas  weight  flow  ratio,  and  is  input  in  the 
TD2  namel ist. 
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SECTION  6 


NEW  OUTPUT  FOR  CFC 


The  output  from  the  Chamber  Flowfield  Code  is  similar  to  the  output 
from  the  FCT  code  except  for  some  new  tables  and  error  messages.  Most  of 
the  new  output  is  self-explanatory.  The  major  new  tables  and  messages  are 
described  below. 

Particle  Impact  Conditions 

Impact  conditions  are  printed  automatically  at  each  output  step  for 
each  particle  size  group  impinging  upon  the  nozzle  walls.  An  example  of  the 
printout  is  shown  in  Table  2.  The  variables  in  the  table  are  described 
separately  in  Table  3.  The  conditions  required  to  calculate  erosion  rates, 
such  as  particle  velocity,  angle  of  impact,  and  impingement  mass  flux,  are 
printed  in  the  output  table  along  with  the  impact  locations.  The  particle 
diameter  is  obtained  from  the  particle  group  size. 

Most  of  the  values  in  the  output  table  are  nondimensional .  To  convert 
them  to  engineering  units  they  must  be  multiplied  by  the  normalizing 
quantities  indicated  in  Table  3.  The  values  of  the  normalizing  quantities  are 
printed  near  the  beginning  of  the  output. 

Particle  Streamline  Details 

Details  of  particle  streamline  relative  positions,  mass  flow 
between  adjacent  streamlines,  and  other  related  variables  can  be  printed 
out  by  the  user,  if  desired.  By  setting  PTPART  to  .True.,  two  additional 
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Table  2.  Example  of  Particle  Impact  Conditions  Table  Printout 
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Table  3.  Description  of  Variables  in  Particle  Impact  Conditions  Table 


Variable 

Description 

Units 

Normalizing 

Quantity 

KK 

Sequence  number 

-- 

-- 

L 

Particle  streamline  designation  (old)* 

-- 

-- 

I 

Mesh  point  counter  along  transformed 
axial  coordinate 

-- 

-- 

X 

Distance  along  the  transformed 
axial  coordinate 

none 

Z 

Distance  along  the  physical  axial 
coordinate 

none 

Throat 

radius  (RSI) 

R 

Distance  along  the  physical  radial 
coordinate 

none 

Throat 

radius  (RSI) 

U 

Particle  impact  velocity  component 
in  the  axial  direction 

none 

AN 

V 

Particle  impact  velocity  component 
in  the  radial  direction 

none 

AN 

VTOT 

Magnitude  of  particle  impact  velocity 
vector 

none 

AN 

THETA  IMP 

Particle  impingement  angle 

degrees 

-- 

ROP 

Particle  spatial  density  at  impact 

1  ocation 

none 

RHON 

FLUXMP 

Particle  impingement  flux 

AN  X  RHON 

*New  and  old  particle  streamline  numbers  are  described  later  in  this 
section . 


tables  are  printed  for  each  particle  size  group  at  each  output  step.  The 
first  table  gives  the  L(old)-L{new)  matrix  shown  in  Table  4.  The  L(old), 

L(new)  particle  streamline  labeling  scheme  is  explained  in  detail  in  Volume  I 

of  the  Particle  Impact  Erosion  Final  Report  (Reference  1).  Essentially,  the 
L(old)  streamline  label  is  applied  sequentially  to  the  streamlines  at  the 
inflow  boundary,  starting  from  1  at  the  axis  to  NIPPG  at  the  nozzle  wall.  The 

L(old)  identifier  remains  fixed  for  each  streamline.  In  comparison,  the 

L(new)  identifier  depends  upon  the  streamlines  local  position.  At  each 
constant  X  mesh  location  (i.e.,  at  each  value  of  I,  the  mesh  counter  in  the  X 
direction),  the  streamlines  are  renumbered  starting  from  1  at  the  axis  and 
increasing  with  Y.  The  L(new)  labels,  therefore,  define  the  relative 
positions  of  the  streamlines  at  a  given  location,  and  these  positions  change 
with  the  occurrence  of  trajectory  crossovers. 

In  cases  where  steamlines  are  terminated,  the  numbering  scheme 
described  above  is  slightly  modified.  Streamline  terminations  occur  when 
subsequent  streamlines  impact  the  nozzle  walls.  In  this  case,  the  terminated 
streamlines  are  assigned  to  the  top  of  the  L(new)  table  and  are  removed  from 
further  calculations.  Streamlines  are  also  terminated  if  they  return  to  the 
chamber  wall,  or  if  the  trajectory  calculations  fail.  In  such  cases,  the 
streamlines  are  assigned  to  the  bottom  of  the  L(new)  table  from  the  outset, 
and  are  removed  entirely  from  the  particle  density  calculations.  However, 
particle  mass  flux  is  conserved  by  redistribution  among  the  remaining 
streamtubes. 

An  example  of  the  second  table  printed,  when  PTPART  is  set  to  true,  is 
shown  in  Table  5.  The  various  quantities  used  in  the  particle  spatial  density 
are  listed.  The  term  VP. DA  is  the  dot  product  of  velocity  and  area  between 
adjacent  streamlines.  DMOOTP  is  the  mass  flow  rate  between  adjacent 
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Table  5.  Example  of  ROP-YP  and  Other  Arrays  Printout 


streamlines,  and  ROP  represents  the  particle  spatial  density  in  that  region. 
The  last  term,  YP,  is  the  location  of  the  particle  streamline  in  the 
transformed  Y  coordinate. 

The  table  also  indicates  the  presence  of  wall  impact  or  trajectory 
crossover,  at  each  value  of  I,  by  the  variables  IMPW  and  ICR.  The  maximum  and 
minimum  values  of  streamlines  are  also  printed,  indicating  the  number  of 
streamlines  terminated.  The  number  of  entries  for  VP. DA  generally  exceeds  the 
maximum  number  of  streamlines.  The  extra  location  is  used  to  store  wall 
impact  conditions. 

Trajectory  Integration  Details 

Details  of  each  iteration  of  the  particle  trajectory  calculations  can 
be  printed  out  by  the  user  by  setting  PTDIRJ  to  .True.  An  example  of  the 
output  is  partially  reproduced  in  Table  6.  The  output  shown  is  for  the  first 
particle  group  (K  =  1)  and  the  first  streamline  (L  =  1).  The  iteration  number 
is  given  in  the  first  column.  The  rest  of  the  variables  have  their  usual 
meanings.  Variables  followed  by  1  refer  to  the  starting  values  of  those 
variables  for  each  set  of  iterations.  Variables  followed  by  2  refer  to  their 
values  after  each  iteration. 

Since  the  total  printout  of  this  table  is  massive,  it  should  be  used 
mainly  for  debugging  purposes.  It  can  also  be  used  to  check  the  details  of 
particle  trajectories,  since  the  regular  printout  lists  streamline  positions 
only  at  each  constant  X  mesh  location.  If  the  particle  crosses  the  same  mesh 
line  twice,  the  particle  conditions  at  those  locations  are  written  over  and 
some  portions  of  the  trajectory  may  go  unrecorded.  This  is  especially  true 
for  particles  which  cross  the  axis,  and  is  the  reason  why  particle  streamline 
plots  sometimes  look  jagged  and  discontinuous.  It  should  be  emphasized  that 
although  intermediate  particle  locations  and  conditions  may  not  be  stored  and 
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Table  6.  Example  of  Particle  Integration  Calculations  Table 


f 
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plotted,  the  particle  trajectories  are  calculated  correctly,  as  may  be 
verified  by  printing  the  detailed  trajectory  table. 

Miscellaneous  New  Tables 

The  code  prints  the  WPSCOL  and  DMDOTP  arrays  at  the  initial  line  for 
each  particle  group  at  each  output  step.  An  example  is  shown  in  Table  7. 

The  variable  WPSCOL  was  described  earlier.  If  the  option  to  calculate  WPSCOL 
internally  by  specifying  WPEXP  is  employed,  the  calculated  values  of  WPSCOL 
will  appear  here.  Since  the  calculated  values  may  vary  slightly  at  each 
iteration  step,  they  are  not  printed  toward  the  front  of  the  output  as  are  the 
input  values  of  WPSCOL.  The  OMOOTP  array  at  the  initial  line  is  different 
from  the  array  described  earlier.  Here  the  variable  refers  to  the  particle 
mass  flow  rate  between  adjacent  mesh  points  at  the  chamber  boundary. 

The  code  also  prints  a  table  toward  the  front  of  the  output  which  lists 
the  inflow  boundary  conditions  at  each  mesh  point  along  the  chamber  boundary. 
An  example  is  given  in  Table  8.  These  may  differ  from  the  input  values, 
since  the  mesh  point  values  are  obtained  by  interpolation  along  the  boundary 
length.  The  table  may  be  used  to  check  the  input  and  interpolated  values  of 
the  boundary  conditions. 

Error  Messages 

The  code  prints  error  messages  if  it  encounters  problems  during 
execution.  Errors  occur  most  commonly  during  particle  trajectory 
calculations,  examples  of  which  are  shown  in  Table  9.  In  CFC,  the  effect  of 
errors  is  generally  not  as  catastrophic  as  in  FCT.  Logic  was  introduced  in 
the  code  to  minimize  the  problems  and  inaccuracies  associated  with  the  errors, 
and  to  continue  with  the  main  body  of  calculations.  The  first  error  message 
in  Table  9  serves  as  an  example  of  this  procedure.  It  indicates  a  failure  to 
converge  in  the  particle  trajectory  integration  calculations.  In  order  to 
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Table  9.  Examples  of  Particle  Trajectory  Calculations  Error  Messages 


increase  the  chances  of  convergence,  the  number  of  iterations  have  been 
increased  to  30  in  CFC.  Also,  it  was  found  that  even  in  cases  where  the 
integration  had  not  formally  converged  in  30  iterations,  the  values  were  very 
close  to  the  final  solution.  Hence,  the  particle  trajectory  calculations  in 
CFC  continue  forward  from  the  last  iteration  step  even  if  the  integration  has 
not  fully  converged.  A  message,  however,  is  printed  to  inform  the  user  of 
this  occurence. 

The  second  error  message  in  Table  9  is  printed  when  a  streamline 
encounters  problems  in  the  XYLOC  subroutine,  which  calculates  the  position  of 
the  streamlines  on  the  X-Y  mesh.  If  the  local  mesh  geometry  is  poor,  the 
streamline  may  become  trapped  within  a  mesh  cell.  In  this  case,  the  entire 
streamline  is  terminated  by  the  code  {cf.,  earlier  discussion).  The  particle 
mass  flow  associated  with  the  streamline  is  distributed  among  adjacent 
streamlines,  so  that  the  total  particle  mass  flow  is  conserved.  Hence,  the 
code  calculations  are  unaffected  by  the  error  except  for  a  degradation  in  the 
particle  flowfield  resolution. 

The  third  error  message  in  Table  9  occurs  when  a  streamline  returns 
to  the  inflow  boundary.  Presently  the  code  cannot  handle  these  situations, 
and  the  streamlines  are  terminated  as  with  XYLOC  errors. 

The  final  error  message  occurs  when  two  or  more  streamlines  impact  the 
nozzle  wall  between  two  adjacent  mesh  points,  a  situation  the  code  logic 
cannot  currently  handle.  If  this  error  occurs  during  the  final  time  step  of  a 
run,  it  may  result  in  problems  since  the  impact  conditions  will  not  be  known. 
The  user  may  attempt  to  eliminate  the  error  by  rerunning  the  case  with 
different  values  of  WPSCDL  for  the  impacting  streamlines. 
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SECTION  7 

SAMPLE  PROBLEM  USING  CFC 

The  Hercules  NTF  001  nozzle,  fired  on  a  Brutus  motor,  was  used  to 
demonstrate  the  subsonic-transonic  total  recession  prediction  methodology 
developed  for  the  Particle  Impact  Erosion  Program.  Details  of  the  firing,  the 
actual  chamber  and  nozzle  geometries,  and  the  approximations  made  for  the  CFC 
calculations  are  given  in  Reference  1.  The  CFC  results  at  an  early  burn  time 
are  presented  here  as  a  sample  problem. 

The  input  data  list  for  the  sample  problem  is  shown  in  Table  10.  The 
properties  previously  communicated  to  FCT  from  ODE  were  obtained  from  ACE 
(Reference  9)  calculations.  A  summary  of  the  ACE  results  is  shown  in 
Table  11.  Average  properties  between  the  chamber  and  throat  were  used  in  CFC. 
Also,  since  the  expansion  in  CFC  extends  only  to  the  transonic  region,  the 
expansion  ratio  (EPM),  gas  velocity  (UGM),  and  particle  temperature  (TPM)  were 
selected  at  a  location  where  the  Mach  number  was  slightly  greater  than  one, 
rather  than  where  the  particle  melt  temperature  occurs.  (See  column  labeled 
location  "M"  in  Table  11).  The  particle  latent  heat  (DHM)  was,  therefore, 
dummied  by  selecting  a  very  small  (negligible)  value.  The  TPS  and  UGS  values 
were  selected  at  a  higher  Mach  number  (Column  labeled  location  "S"  in 
Table  11). 

Selected  outputs  are  shown  in  Table  12.  The  particle  impact 
conditions  and  other  parameters  are  shown  at  the  1,000th  time  step.  The 
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Table  10.  CFC  Sample  Problem  Input 


TITLE  TEST  CASE  lA  --  NTF-001  /f«0  BURN  T ME  6E0N 

SPP 
tSPP 
TD^Pil.  , 

RSI  =  ^.2*10. 

INLET  =  3,  PHI  =  RR.Ot  RUT  U::! .  350« 

^I  =  -l.Slllt  RI  =  2.0533, 

IUALL=3,  RUT0=1.35a,  THeTA=lA.40,  ePS=3.a, 

CNNT  =  0.1,  0.3,  0.5,  0.7*  0.9*  1.0,  1.1,  NCHNT  =  7, 

CPRT  =  1.6,  1.5,  1.4,  1.3,  1.1,  1.0,  O.R,  0.6,  0.4,  NCPRT=9, 
PLTT=1,  PRUIPL=1,  NPLTS=2»1*0,3*1 ,  NSPL T 3=2 • 1 , 0 , 3 *1 , 
tENO 

•ro2 

PR  -  0.41125,  6NS0  =  0.46646E-04,  CAPR  =  0*^607, 

CPS  =  R49d.0,  CPL  =  8498.0, 

TPH  =  4759.1,  OMN  =  l.OE-06,  EPH  =  1.1600,  UGN  =  4325., 

TPS  =  4573.9,  UGS  =  4808., 

PC  =  S3S.0,  TOO  =  5497.8* 

UPU6T  =  0.50693, 

REG2  3  .TRUE.,  PTRAjIP  s  «TRUE.t  /NAXArSfOPt 
NX  =  31,  NT  =  21, 

PKK  =  105*0.2,  PLK  =  105*0.5, 

ITHTN=1,  POIN=21*1.0,  TOIN=?l*1.0, 

THIN  2  i.0C-Q6*19*-89.0*-92.0, 

IPCALCsl* 

THPK(1,1|s  0.0,19*>a9.Of92.0* 

THPK(1*2I=  0.0,19*-a9.a,-92.0, 

THPK<1,3»2  0.0,19**89.0, -92.0, 

NPRT23*10a0,  NPLT23*1000, 

OENO 

aT03 

ROOT  s  0.3557, 

ALFAAP  2  0.70, 

IPATCH  2  3, 

2PATCH  2  -8.0000,-4.4056,-0.81127, 

RPATCH  =  0.000*  0.5997*  1.19931* 

ICHNB211*  ICHN8228,  ICHMBls4, 

ZCHNB  2  -47,619, -47*6ia*-4T*6t7,-47. 616, 

-4  3.172,-38.728,-34.2  84,-29.8  40, 

-19.170,-10.340,  -1.5107, 

RCH8M  2  J.JOO,  1.038,  2.076,  3.114, 

3.113,  3.112,  3.111,  3.110, 

2.765,  2.420,  2.0743, 

FPERCH  2  11*1  .0  , 

BEXPCH  2  11*0.343,  C8HNCM  2  11. 0.3874, 

SPCUPL  2,  TRUE  .  ,  K3f’AR22, 

PTPART  2  .TRUF.,  Ni.IIMl  2  ?^,o,  ylHOX  2  0.5292, 

UPEXP  2  l.n,  NPSSUH  25, 

OAHPN  2  .TRUE.,  SHLOI  2  0. 3,  SHLOJ  2  u. 3, 

SEND 
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Table  11.  ACE  Run  Results  Used  in  Sample  Problem 


Location 

Chamber 

Throat 

M 

S 

Pressure  (atm) 

60.45a 

34.615 

20.0 

15.0 

Mach  number 

1.0 

1.433 

1.624 

Temperature  ("R) 

5.497.8b 

5.120.4 

4,759.1 

4,573.9 

Particle  enthalpy  /  Btu  \ 

-5,045.6 

-5,173.8 

Particle  to  gas  mass  ratio 

0.49916 

0.50693 

Gas  viscosity  /  Ibm  \ 
\7t7secj 

0.4665  X  10-4 

0.4450  X  10-4 

Gas  Prandtl  number 

0.41042 

0.41209 

Velocity  /ft  \ 

[sec/ 

3,127 

4,325 

4,808 

Expansion  ratio 

1.000 

1.160 

^Average  chamber  pressure  =  888.6  psia 

^Flame  temperature  (heat  loss  to  exposed  insulation  in  aft  closure  considered) 
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Table  12.  (Continued) 


Table  12.  (Concluded) 
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calculations  had  not  entirely  converged  by  that  time,  still  having  small, 
residual  transients.  However,  it  has  been  found  that  the  flowfield  closely 
approaches  steady-state  in  approximately  1,000  steps.  Hence,  a  limit  of 
1,000  steps  is  usually  sufficient  for  calculating  particle  impact  conditions. 
The  slow  rate  of  final  convergence  is  almost  certainly  caused  by  the  low  Mach 
numbers  near  the  chamber  boundaries.  Substantial  improvements  in  the  rate  of 
convergence  would  require  the  development  of  new  calculation  procedures  in  the 
low  Mach  number  regime. 

Figure  5  shows  the  CFC  mesh  geometry  for  the  sample  problem,  and 
Figure  6  shows  the  trajectories  for  the  three  particle  size  groups  evaluated. 
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Figure  5.  CFC  Mesh  Geometry  for  the  NTF  001  at  0  Sec  Burn  Time 
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(a)  Particle  Size  Group  1 
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(b)  Particle  Size  Group  2 
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(c)  Particle  Size  Group  3 


Figure  6.  Particle  Trajectories  for  the  NTF  001  Demonstration 
«  '  at  0  Sec  Burn  Time 
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SECTION  8 


CMAE  ANALYTICAL  BASIS 

Particle  erosion  is  coupled  with  gas-surface  thermochemical  ablation  in 
the  no22le  entrance  and  throat  regions.  A  comprehensive  procedure  already 
existed  to  predict  recession  in  rocket  nozzles  due  to  ablation  alone 
(References  12  and  13).  The  procedure  was  modified  to  account  for  the  effects 
of  impact  erosion.  It  was  assumed  that  particulate  erosion  did  not  affect  the 
thermochemistry,  kinetics,  and  boundary  layer  aspects  of  the  calculation 
procedure.  Hence,  the  onlv  modifications  required  were  to  the  recession 
calculations.  The  CFC  provides  the  impact  parameters  which  are  input  to  the 
CMAE  code. 

The  base  code  for  the  recession  calculations  was  the  Charring  Material 
Ablation  (CMA)  code  (Reference  5).  The  new  code  calculates  mass  loss  due  to 
both  erosion  and  ablation;  hence,  it  was  called  the  Charring  Material 
Ablation-Erosion  (CMAE)  code.  The  major  modifications  to  the  code  were 
inclusion  of  a  subroutine  to  calculate  the  erosion  contribution,  and 
modification  to  the  surface  energy  balance  to  include  erosion  effects.  The 
erosion  contribution  is  calculated  using  the  model  for  mechanical  and  chemical 
erosion  described  in  Reference  1.  The  modification  to  the  surface  energy 
balance  is  discussed  below.  Again,  It  is  assumed  that  the  reader  is  familiar 
with  the  basic  concepts  and  procedure  (References  12  and  13). 
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Erosion  Model 

The  erosion  model  in  CMAE  is  in  the  form  of  a  G-law  with  two 
components:  mechanical  erosion  due  to  mechanical  mass  removal,  and  chemical 
erosion  due  to  chemical  reactions  between  molten  particles  and  the  surface. 

G,  which  is  the  ratio  of  mass  flux  of  material  removed  to  mass  flux  of 
incident  particles,  can  be  related  to  particle  impact  conditions 
(Reference  11),  as  follows: 

G  =  a  (sin  a)"*  T®  (22) 


where 


Vp  =  particle  impact  velocity 

Dp  =  diameter  of  impacting  particles 

a  =  angle  of  impact  relative  to  the  surface 

T  =  target  temperature 

a,  b,  c,  d,  e  are  empirically  determined  constants 

Each  of  the  mechanical  and  chemical  components  of  erosion,  G^gj-h  and 
Gchem  3''®  ’f’  Equation  (22);  the  total  erosion  being  given  by  the 

sum  of  the  two  components.  A  special  form  of  the  temperature  dependence  was 
formulated  for  rocket  nozzle  erosion.  The  limited  data  available  on 
temperature  indicated  that  no  erosion  occurred  below  3,000°R  (T^)  and  that 
erosion  was  essentially  independent  of  temperature  above  4,500°R  (T2).  The 
functional  dependence  of  erosion  on  temperature  between  these  two  limits  was 
assumed  to  be  smooth  and  was  modeled  by  the  hyperbolic  tangent  function: 
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(23) 


where 


T 

0 


Tj  +  Tp 

-  ,  the  midpoint  between  the  temperature  limits 


AT  =  n(T2  -  Tj),  a  stretching  parameter.  (The  recommended 
value  of  n  is  1/4.) 

The  G-law  used  in  the  CMAE  code  incorporates  both  the  power  law  and 
the  hyperbolic  tangent  forms  of  temperature  dependence,  as  follows; 

G  =  a  (sin  a)"  f  f(T)  (24) 

The  user,  therefore,  has  the  option  of  selecting  the  hyperbolic-tangent 
function  (by  setting  e  =  0) ,  or  the  power-law  function  (by  setting  Tg  =  0), 
or  any  combination  of  the  two. 

Since  the  combined  effect  of  erosion  and  ablation  is  important,  the  CMA 

surface  energy  balance  was  modified  to  account  for  the  coupling  between  the  two 

mechanisms.  As  shown  in  the  following  section,  he  surface  energy  balance  is 

modified  only  by  one  term:  L  irip  AHn,  where  liip  is  the  mass  flux  removed 

®cnem  ^  ®chem 

due  to  chemical  erosion,  AHr  is  the  heat  of  reaction  and  E  signifies  summation 

over  all  particle  types  and  sizes.  The  value  of  irip  is  calculated  from  the 

chem 

G-law  for  chemical  erosion,  and  aHr  is  an  input  variable. 

The  addition  of  this  new  term  affects  the  other  energy  balance 
parameters.  In  particular,  the  conduction  term  is  directly  affected  since  it 
is  influenced  by  the  recession  rate.  It  is,  therefore,  important  to  account 
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for  the  couple  effects  of  ablation,  erosion,  and  in-depth  conduction  in 
calculating  the  total  recession  in  the  entrance  and  throat  regions  of  rocket 
nozzles. 

Surface  Energy  Balance 

The  surface  energy  balance  (SEB)  was  modified  in  the  new  (CMAE)  code  to 
include  the  energy  fluxes  associated  with  particle  erosion  effects.  The 
formulation  is  similar  to  the  one  used  in  the  ABRES  Shape  Change  Code  (ASCC) 
(Reference  11).  The  following  sketch  illustrates  an  ablating-eroding  surface, 
receding  with  the  rate  s  in  the  surface  normal  direction.  A  differential 
control  volume  is  fixed  to  the  surface.  The  important  energy  fluxes, 
including  erosion-related  terms,  to  and  from  the  surface  are  shown. 


^diff  **rad  **rad  ^^chem^ 


in  out 

1  1  1 

1  Particle  mass 

---l-.i _ L 

- H - 

loss  energy 

Nozzle 

surface 


"cond  ^3",  <’tc  * 


The  energy  inputs  to  the  control  volume  include  the  following: 

•  The  diffusive  energy  fluxes  from  the  boundary  layer,  q^iff 

•  The  radiative  energy  absorbed  by  the  surface, 

in 

•  The  energy  released  due  to  chemical  interactions  of  the  particles 

with  the  surface,  l  nig  AHp,  where  AHr  is  the  heat  of  reaction 

chem 

(ilg  is  the  mass  flux  removed  due  to  chemical  erosion,  and  the  Z 
chem 

signifies  summation  over  all  particle  types  and  sizes 


f 


« 


t 


1 

•  The  enthalpy  flux  carried  to  the  surface  by  the  pyrolysis  gases, 

trig  hg,  where  trig  and  hg  are  the  mass  flux  and  enthalpy,  respectively, 
of  the  pyrolysis  gases 

•  The  energy  flux  due  to  the  relative  motion  of  the  control  volume, 

(m^c  '"e)  ^s»  where  is  the  mass  flux  of  the  surface  removed 
due  to  thermochemical  ablation,  mg  is  the  mass  removed  due  to 
erosion,  and  hg  is  the  sensible  enthalpy  of  the  solid  material  at 
the  surface 

The  energy  fluxes  leaving  the  control  volume  include  the  following: 

•  The  energy  flux  carried  into  the  material  due  to  conduction, 

‘Icond 

•  The  radiative  energy  flux  emitted  by  the  surface, 

out 

•  The  energy  carried  away  due  to  gross  blowing  at  the  surface,  (pv)y,h„, 
where  (pv)^,  is  the  blowing  mass  flux,  (pv)y^  =  m+c  +  itiq,  and  h„  is  the 

^  i 

enthalpy  of  the  blown  gases  at  the  surface  | 

•  The  energy  flux  lost  from  the  surface  due  to  erosive  mass  removal 

The  last  item  equals  the  product  of  the  mass  flux  removed  due  to  erosion, 

mg,  and  the  sensible  enthalpy  of  the  surface  material,  hg.  Hence,  the  surface 
energy  balance  can  be  written  as: 

"diff  *  ■’rad  ‘"r  *  "gS  * 

in 

=  +  (iTi*  +  m  )h  +  m  h  (25) 

^cond  ^rad  '  tc  g'  w  e  s  '  ' 

out 

The  final  terms  on  each  side  of  the  equation  cancel  out,  leaving  the  chemical 
interactions  term  as  the  only  new  term  in  the  modified  surface  energy  balance. 
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However,  erosion  influences  most  of  the  remaining  terms  in  the  SEB. 

Conduction  is  directly  affected,  since  conduction  in  a  moving  medium  depends 
upon  the  rate  of  movement  of  the  medium.  The  other  terms  are  indirectly 
affected  since  they  must  be  adjusted  in  order  to  maintain  an  energy  balance  at 
the  surface. 

The  modified  SEB  solves  for  the  wall  temperature  and  the  recession 
rate.  An  interative  procedure  is  used  to  solve  implicitly  for  the  unknowns  at 
each  time  step.  The  procedure  is  identical  to  the  one  used  in  the  base  (CMA) 
code  except  that  the  erosion  calculation  routine  is  included  in  the  implicit 
procedure. 
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SECTION  9 


\ 


CMAE  INPUT 


This  section  describes  only  the  new  input  required  for  CMAE,  which  is 
inserted  between  the  CMA  surface  time-dependent  boundary  conditions  and  the 
surface  thermochemistry  tables.  The  erosion  input  contains  a  control  card 
describing  the  number  of  impacting  particle  groups,  the  number  of  particle 
types,  and  the  names  of  those  particle  types.  If  no  erosion  is  desired,  a 
blank  card  must  be  input  to  denote  zero  particles.  After  the  control  card, 
the  following  data  are  input:  a  set  of  particle  time  histories  for  each 
particle;  the  coefficients  for  the  mechanical  erosion-mass  loss  law;  and  the 
coefficients  for  the  chemical  erosion-mass  loss  law  and  the  temperature 
dependent  heats  of  reaction  for  each  particle  type.  The  formats  for  each  of 


the  card 

set  types 

are  described  below: 

Control 

Ca  rd 

Col umn 

Format 

Variable  Description 

Units 

1-5 

15 

Number  of  particle  groups  (max  =  10) 

None 

6-10 

15 

Number  of  particle  types  (max  =  2) 

None 

11-15 

A5 

Alphanumeric  description  of  particle  type  1 

None 

16-20 

A5 

Alphanumeric  description  of  particle  type  2 

None 
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Particle  Group  Time  History 

One  card  set  for  each  particle  group  is  required. 

Card 


Number 

Col umn 

Format 

Variable  Description 

Units 

1 

1-5 

15 

Particle  type  (1  or  2) 

None 

11-20 

FIO.O 

Diameter  of  particle 

ym 

2+ 

1 

11 

Flag,  nominally  blank,  punched  to 
indicate  last  card  of  time  table 

None 

2-10 

F9.0 

Time 

sec 

11-20 

FIO.O 

Impacting  mass  flux 

Ibm/ft^-sec 

21-30 

FIO.O 

Particle  impact  velocity 

ft/sec 

31-40 

FIO.O 

Particle  impact  angle 

deg 

Mechanical  Erosion-Mass  Removal  Law 

These  coefficients  correspond  to  those  in  Equations  (22)  and  (23)  for 
mechanical  erosion. 


Col umn 

Format 

Variable  Description 

Units 

1-10 

FIO.O 

a 

None 

11-20 

FIO.O 

b 

None 

21-30 

FIO.O 

c 

None 

31-40 

FIO.O 

d 

None 

41-50 

FIO.O 

e 

None 

51-60 

FIO.O 

To 

°R 

61-70 

FIO.O 

at 

°R 

Chemical 

Erosion-Mass 

Removal  Law  and  Energy 

One  card  set  for  each  particle  type  is  required.  The  coefficients 


correspond  to  those  in  Equation  (23)  for  chemical  erosion.  (The  same 
temperatures,  Tg  and  aT,  are  used  from  above.) 
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Col umn 

Format 

Variable  Description 

Units 

1-10 

FIO.O 

a 

None 

11-20 

FIO.O 

b 

None 

21-30 

FIO.O 

c 

None 

31-40 

FIO.O 

d 

None 

41-50 

FIO.O 

e 

None 

1 

11 

Flag,  nominally  blank,  punched  to 
indicate  last  card  of  temperature 
table 

None 

2-10 

F9.0 

Temperature 

“R 

11-21 

FIO.O 

Heat  of  reaction  for  the  following 
reaction: 

Btu/lb-carbon 

n  Particle  +  C* — ►Product  +  aHr 

SECTION  10 


SAMPLE  PROBLEM  USING  CMAE 


The  calculation  for  total  recession  in  the  nozzle  entrance  region  of 
the  NTF  001  firing  is  presented  here  as  a  sample  problem.  The  firing  results 
were  used  to  demonstrate  the  subsonic-transonic  total  recession  prediction 
methodology  developed  for  the  Particle  Impact  Erosion  program,  as  described  in 
Reference  1. 

Table  13  presents  the  CMAE  input,  while  Table  14  illustrates  the  new 
output  of  erosion  data  and  recession-temperature  at  selected  times. 


Table  13.  CMAE  Sample  Problem  Input 


CMA  ANALYSIS  OF 

PITCH  FIBER  1 

FORMARD  LOCATION  THROAT 

THERMOCHEMISTRY 

BASED  ON  HTPB 

A  ' 

1. 

0. 

B 

1. 

0. 

C 

220. 76 

120.76 

to 

26 - 0.0 

63.6 

1.0  25.0 

.006 

2 

530. 

.002 

Z 

- 530i - 

.006 

2 

530. 

.006 

2 

530. 

.008 

■  Z' 

'330. 

.010 

2 

530. 

.012 

2 

530. 

.016 

-2- 

- rsoi - 

.016 

2 

530. 

.018 

2 

530. 

.020 

-  -2- 

- 530. 

.030 

2 

530. 

.060 

2 

530. 

.050 

-2 

590s - 

-  -  -iOOO 

2 

530. 

.  170 

2 

530. 

.20 

2 

330o 

2 

530. 

.25 

2 

530. 

.25 

~  390« - 

■*  ■ 

3 

530. 

.25 

3 

530. 

.25 

3 

330« 

3 

530. 

.25 

530i - it» 

.01339 

960.  .31 

.01105 

1660.  .39 

1360.- 

.008969 
- .007692 

2,25  IN  RADIUS 


1.68  IN  THICK 
2  MIL  ROUGH  NTTOOIER 


1. 

1. 

1. 

.1-  - 

0. 

2.91518 


1. 

1. 

1. 

0. 


1. 

1. 

1. 

2.5 

0. 


90000. 

90000. 

90000. 

0. 


2660. 

2960. 


.665 

.685 


.006J62 

.005556 

.005161- 


kfllS 

.1115 

.ni5 

-.?115 

.T115 

.7115 

-.M19 


3960. 

.515 

.005069 

.7115 

6660. 

.520 

.005067 

.7115 

- - 

— i926  — 

.005067 

.7115 

-1  7660. 

.526 

.005067 

.7115 

530. 

.15 

.01339 

.7115 

--  -960. - 

-.31 

.01105 

.7115 

1660. 

.39 

.008969 

.7115 

1960. 

.63 

.007697 

.7115 

-  --  2660. - 

.665 

.006362 

.7115 

2960. 

.685 

.005556 

.7115 

3660. 

.505 

.005161 

.7115 

3960. 

.515 

.005069 

.7115 

6660. 

.520 

.005067 

.7115 

6960. 

.526 

.005067 

.7115 

--t-T660^. 

-.526  - 

.005067 

.7115 

.6 

536. 
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Table  13.  (Continued) 


— T-I08.-Z~ 

- -  --  - 

500. 

.210 

.760-4 

1000. 

.362 

.750-4 

- 1300i  - 

- .959 

.380-4 

2000. 

.735 

.650-4 

3000. 

.687 

.160-4 

~4000. 

.630 

.390*3 

5000. 

.590 

.735-3 

6000. 

.590 

.136-2 

• 

"  j590 - 

•  i 

l  0.0 
0.0 

- 

10.0 

1  63.4 

- 3 - 

1 

0.0 

— 2<r.t> 

40.0 

1  65.0 
- 1 - 

0.0 

20.0 

I  65^0 
1 

- 0.0 

20.0 
40.0 
X  65.0 


7000. 

0.0 

-t43c^rf— 

-1436.27 

-1436.27 

-t - 6t - 

90.12 
9.91 

- a.-oo 


-4tO-.-94- 

416.54 

416.54 


615.87 


-ti-50925- 

1.56425 

1.56425 


55.569 

55.569 

55.569 


.641E 

.13356 

-900. 

1650. 

1655. 

4970. 


0.00 

0.09 

12.01 

3.67 

0.00 

204.85 

tt.96— 

14.87 

0.00 

0.00 

.7987 

.8 

-24500. 

-24200. 

”Z9OP0« 


615.87 

615.87 


568.54 

541.12 

541.12 


52tv4t- 

527.07 

527.07 

1.236 

.8 


22.174 

tr.  if4— 

22.174 

22.174 

23.6805 

12.582 

12.582 

»i65a 
12.825 
12.825 
h.029- 
.2992 
.8 


3750. 


375. 


4975.  -12500. 

1  7200.  -12250. 

— -i68499 - trO - 

55.5690  0.0  0.0 

55.5690  0.0  0.0 


1 


500. 

1000. 


55.5690  0. 
55.5690  0. 
55.5690  0. 
55.5690  0. 
55.5690  0. 

.569- 
55.5690  1. 
55.5690  1. 
55.5690  1. 
55.5690  1. 
55.5690  1. 


0. 

0. 

-0i 

0. 

0. 


07056  0 
.07056  0 
.07056  0. 
.07056  0. 
.07056  0 
.Ot056-«. 


2000. 
2500. 
5000. 
3500. 
4000. 
298. 
00002  500. 
000031000. 
000421200 
001261400 
002481600 
004091800. 


00000.667  -825.176-1899.741-1 
00000.667  -561.750-1726.032-1 
00000.^69—270.856-1537.637-1 
00000.667  '  20.489-1338.908-1 
00000.667  335.680-1037.704-1 
00000.667 — 662.455— 822.667-1 
00000.667  995.749  -604.988-1 
00000.667  1334.797  -385.024-1 

.667—1800. - —1800.  1 

00000.667-1559.878-1467.181  1 
00000.667-1075.703-1114.499  I 
80000.667  -693.430  -819.564  1 
00000.667  -472.767  -645.031  1 
00000.667  -J40.633  -540.906  1 
00000.667  — 212.769  -442.036  I 


CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

CHAR 

C4 

C6 

C* 

C4 


0.5010400 
0.5010400 
0.5010400 
0.5010400 
0.5010400 
0.5010400 
0.5010400 
0.5010  00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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Table  13.  (Concluded) 


15.5670  1.0709S~ir.a06aa2000.00000.66T  — ^TtTin--=3«.^62  1 

55.5690  1.07056  0.000532200.00000.667  56.993  -237.973  1 

55.5690  1.07056  0.012102600.00000.667  206.662  -120.730  1 

-51.^690-nv0709r  0.019602600.00000.667  377.227-  --6^020  1 

55.5690  1.07056  0.060792000.00000.667  612.200  150.910  1 

55.5690  1.07056  0.103963000.00000.667  1003.606  626.667  1 

-53.5690  1.0T056  0.160303200.00000.667  1372^003— 730i 735  1 

55.5690  1.07056  0.233213600.00000.667  1007.367  1123.301  1 

55.5690  1.07056  0.363623600.00000.667  2303.177  1671.913  1 
-35.5630  1.07050  0.6056130e0«00000^667-3093^<09  260U700- 1 

55.5690  1.07056  1.013006000.00000.667  3020.351  3192.795  1 
55.569  0.00000  0.00001  290.  .667-1000.  -1000.  1 

55.5690  0.00000  0.00002  500. 

55.5690  0.00000  0.000031000. 

55.5690  0.00000  0.000621200. 

55I5690  OIOOOOO  o!o0260160o! 

55.5690  0.00000  0.006091000. 

55.5690  0.00000  0^006002000; 

55.5690  0.00000  0.000532200.00000.667  56.393  -237.973  1 

55.5690  0.00000  0.012102600.00000.667  206.662  -120.730  1 


00000.667-1075.703-1116.699  1 
00000.667  -693.630  -819.566  1 
00000.667  -677.767  -695.031  7- 
00000.667  -360.633  -560.906  1 
00000.667  -212.769  -662.036  1 


C6 

C6 

C6 

C* 

C6 

C6 

C6 

C6 

C6 

C6 

C6 

CHAR 

CHAR 

CHAR 

C6 

C6 

C6 

C6 

-€6  - 

C6 

C6 

C6 


55.5690  0.00000  0.060792000.00000.667  612.280  158.910  1  C* 

55.5690  0.00000  0.103963000.00000.667  1003.606  626.667  1  C6 

35.5690  0.00000  0.180303200.00000*667  1572.009 — 738.733  1  C6 - 

55.5690  0.00000  0.233213600.00000.667  1807.367  1123;301  1  C6 

55.5690  0.00000  0.36362360(^00000.667  2383.177  1671.913  1  C6 

95.5890  0.00000  0.0fl561380^0000»8667  3095.806  2601.100  7  <6 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0  — 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
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Table  14.  CMAE  Sample  Problem  Output 


- EROSION  MODEL  INPUT - 


-pinrrtciT  HO.  t 

TIME 


PARTICLE  TTME  “At-  eAUTtCLE  OIAN.  « 
NASS  FLUX  VELOCITY  ANGLE 


90. 


0.0  9.91000  AlS.aiOOOO  22.174 

-  — - 20.000  -  .0  919.070000  --22i  174 

40.000  -  .0  615.070000  22.174 

69.000  .0  619.070000  22.174 


PARTICLE  NO.  2  PARTICLE  TYPE  AL  PARTICLE  OIAN.  >  131., 

- TIME  NAS9  FLU* - VELOCITY - RNOLE 

U0M/FT2-SECI  IFT/SECI  lOECI 


- t2s9teO - 990i  940000 — 29*900 - 

20.000  3.67000  941.120000  12.502 

40.000  .0  941.120000  12.502 

09.000 - iO - 541.1 20000  -12*902- . 


PARTICLE  NO.  3  PARTICLE  TYPE  AL  PARTICLE  OIAN.  ■  204. 


TINE  NASS  FLUX  VELOCITY  ANGLE 

(LBN/Ff2-SECf  IFT/SECI  I0E6I 


0.0  12.9600  521.410000  .25.693 

20.000  14.0700  527.070000  12.025 

■40*000 - *0  — - 521*070000 — 12*025 

09.000  .0  527.070000  12.829 


NICRONS 


NICRONS 


NICRONS 
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Table  14.  (Continued) 


- G-LAM  CONSTANTS - 

♦i>FAitr  «•  a  •  vELCCiTv  ••  e  ♦  siN«Attetc«  o  ♦  ten^  ••  e 
NOOEL  A  a  C  0  E 
erosion 


- — ^«teaoo-aa  iTaaraa  -liEsaaa - Ttaasee  .o 

CHEMICAL 

.1335000-03  .800000  .800000  .800000  .0 

rnpsrn  «  arso.oo  tnmstm  »  375.00 


TYPE 


AL 


PART  TEMP  ENERCV 

-  -  -  T0E6RI  — tatu/iaN  CA»ai - 


- IIW.WW  - 

1650.00  -26500.0 
1655.00  -26200.0 
-  6670.00  -25000*0- - 

6975.00  -12500.0 
7200.00  -12250.0 


-SURFACE  EOUIIIORIUN  DATA - 


RATIO  OF  NASS  TO  HEAT  TRANSFER  COEFFICIENTS  ■  0.686 

UNEOWAt  OiPPOOICW  ERPONENT  «  0*667 - 

NOMINAL  SURFACE  VIEM  FACTOR  «  1.000  fOPTION  II 
FISSURE  NOOELNCT  USED  FOR  GAS  TERNS 

NO  CHAR  SMELL  CCRRECTION  ON  SURFACE  RECESSION 
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Table  14.  (Continued) 


-Tiwe  sour  — M«o«  surface 

SIEP  IfER  OPfN  RAO  I IRI 
1723  3  1  3.3145 


20.0000  SECONDS  -  -  - 

H  NALL  -  It-COGE  -  HEAfCOEfF  CH/CHO 

(BTU/LBI  IBTU/LB)  ILB/SQ  FT-SECl 
•2AB5.0S  -143*. 27  1.5631  0.99429 


B  li'RlNE  B  PRIME  G 
0.002*7  0.0 


- ABLATION  RATES - 

M  001  CHAR  N  DOT  CAS  N  CHAR  N  GAS 

tLS/S«  ^T-StCI - ItB/ORlG  SO  PTl 

0.002S55  0.0  4.3*3*27  0.0 


SURFACE 


tlNI  /  llN/SECl 
CHAR  10.021 
I. 079 7929/-. 001 15*4 


P7R0LVS1S  I0.9S1 
797929/- .00 1 1 5*4 


RATE 

TOTAL 


- SURFACE  ENER67  FLUX  TERNS - 

- CURRENT  RATES  lBTU/50  FT  SURFACE-Stei - 

ANO  INTEGRATED  VALUES  IBTU/ORIG  SO  FTI 
CONVECTEO  RAOIATEO  RADIATED  CHEMICAL 

- IN - 9N - OUT - OENERATiON 

0.1950*04  0.29*0*03  0.2530*02  0.2340*03 

0.4300*05  0.6310*04  0.4590*03  0.5210*04 


CONDUCTION 

AHAV - 

0.4380*03 

0.9S30*04 


- INTERIOR  ENERGY  TERNS - 

CURRENT  RATES  IBTU/SQ  FT  SURFACE-SECI 

- RNO-INTEORATEO  VALUE*  IBTU/ORIG  SB  9T* 

FVRQL  GAS  OECONF  CONVECTION  STORAGE 

PICK  UP  ABSORPTION  MITN  SOLIDS  IN  SOLID 

—0*0 - - 0.1070*03 - Bt24E0*03^ 

0.0  0.0  0.3*80*04  0.6070*04 


LOSS  AT 
REAR  FACE 

-  0.64*0*01 - 

0.7080*02 


—  TENP—  OENSIfY  ENfHALPT- 
lOEC  Rl  ILB/CU  FTI  IBTU/LBI 
2932.67  120.740  925.79 


- TEMP  -DENSITY  ENTHALOV 

10E6  Rl  ILB/CU  FTI  IBIU/LBI 
2177.72  120.740  574.95 

-2029.19 - 129.740  -  59*»«9 — 


3 

2- 

20BO.A0 

120. 740 

901.0* 

15 

2 

1797.11 

120.740 

410.73 

4 

2 

2B44.B9 

120. 740 

B04.09 

1* 

2 

1558.13 

120.740 

312.75 

- »— 

-2- 

-2891.95 

— 120.740 

0*3.70 

— IP- 

1372.79 

120.740 

242.00  - 

* 

2 

2752.17 

120. 740 

840.05 

18 

2 

12*2.53 

120.740 

203.41 

T 

2 

2*96.02 

120.740 

B13.3B 

19 

2 

1236.15 

120.740 

194.17 

- "B — 

-t- 

—2934.23 

-120.740 

784.03 

—21— 

^10.29 

10B.200 

0.0 

9 

2 

25*7.68 

120.740 

752.42 

22 

3 

533.4* 

108.200 

0.0 

10 

2 

2497.30 

120.740 

718.99 

23 

3 

530.11 

108.200 

0.0 

12 

-r 

2 

-2412.18 

2301.92 

120.740 

120.740 

*79.42- 

*30.53 

24  - 

-9 

-  530.00 

108*200 

0.0 
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Table  14.  (Continued) 


STEP 

1723 

SUftP 

ITER 

3 

PROS 

OPTN 

1 

-  -  20.0000 

SURFACE  H  MALL 

RAO  (INI  (OTu/tai 

3. 3193  •2A03.i>3 

SECONDS 
n  EOOE 

laru/LBi 

•1436.27 

HEAT  COEFF 
Ua/SO  FT-SECI 
1.3631 

CH/CHO 

0.99929 

- CHEMICAL  EAOSION  ENERCV - 

- #41#  i -  - 

-0.201M>»0« 


—NASS 

LOSS  SUNNARV - 

- NASS  Fttt* - 

(LR/SQ  FT-SECI 

*^IIN/SCCl"*^* 

- RECESSION 

UNI 

1 MCHMULVICM ILAL 

ABLATION  0.28330-02 

0.20370-03 

0.46690-02 

*  AM  1 1 LL  E 

NECHANICAL 

EROSION  0.13200»00 

0.13200-01 

0.2271D»00 

PARIICLE 

CHENICAL 

EROSION - »iO301B~dl - 

- 0*04490  02  - 

- 0.1 7200^00 - 

Table  14.  (Continued) 


HEAf  COEFF  CH/CHO 

ILB/SQ  FT-SECI 

1.&4S3  o.«a7a« 


- ABLAflCN  RATES - 

a  FRINE  a  PRIME  C  N  DOT  CHAR  N  DOT  £AS  N  CHAR  N  GAS 

- —  tt«/S(r^T-S«CT - Ita/tHtlG  SG  FT» 

0.04499  0.0  0.047503  0.0  10.873561  0.0 


- - - -  63.4000  SECONDS  -  - 

TWr  WRF  PROB  SURFACE  H  VALt.  -  H  EDGE  - 

STEP  ITER  OPTN  RAO  fINI  (BTU/LBI  IBTU/LBI 

3716  3  1  3.0475  •160B.5S  -1436.27 


SURFACE 

0.9317019/0.0047291 


(INI  /  IlN/SECI 
CHAR  (0.02) 
t.6f960B2/-vOOOOOia  - 


PVROLVSIS  10.98) 

1 .67 96082/-.00000ia 


RATE 

TOTAL 


CONVECTEO 

- IN - 

0.2660*03 

0.9790*05 


- SURFACE  ENERGY  FLUX  TERNS - 

ANO  INTEGRATED  VALUES  (BTU/ORIG  SO  FT) 
RADIATED  RADIATED  CHEMICAL 

-  IN - out - generation 


0.2960*03 

0.2260*05 


0.2200*03 

0.7690*04 


-0.1240*03 

0.0900*04 


PVROL  GAS 
PICK  UP 

-RRTE - OsO - 


TOTAL  0.0 


- INTERIOR  ENERGY  TERNS— 

CURRENT  RATES  lOTU/SQ  FT  SURFACE-SEC) 
lNTeORAfE»^»6tUEfr  ♦B7U/0R1G-SN  ET7 
OECOMP  CONVECTION  STORAGE 

ABSORPTION  WITH  SOLIDS  IN  SOLID 

0*0  0.0790*04  0.1910*05 


CONDUCTION 

AHAV 

0.2110*03 

0.2070*05 


LOSS  AT 
REAR  FACE 

o.iaoo*e2- 

0.7790*03 


NOee  MAT - tOWP  BENS  H Y  «rfHALR» — NOOE  NAT - TENP 

lOEG  R)  ILB/CU  FT)  (BTU/LB)  lOEG  R) 

1  2  5077.73  120.740  2023.45  12  2  4561.76 


3 

4 

-9- 

6 

7 

-B- 

9 

10 


5046.55 

5022.47 

-4971571- 

4954.45 

4910.00 

-4007527- 

4005.90 

4745.15 


120.740 
120.740 
120.740 
120. 740 
120.740 
120. 740 
120. 740 
120.740 


2007.05 

1994.30 

1970520- 

1950.62 

1935.83 

-1909*00- 

1000.93 

1849.15 


14 

15 


17 

21 

-22- 

23 

24 


2  4261.00 

2  3901.12 

. 3705.56 

2  3627.10 

3  1143.76 
-006.90 

538.20 
530.78 


DENSITY 
ILO/CU  FT) 
120.740 
—  120.740 
120.740 
120.740 
120.740 
120.740 
108.200 
108.200 
100.200 
108.200 


ENTHALPV- 
IBTU/LO) 
1753.24 
-  1609.7^  - 
1597.45 
1452.20 
1311.51- 
1271.53 
0.0 

0.0 - 

0.0 

0.0 


Tiir  120.740  -  1000*1 
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Table  14.  (Concluded) 


nur-sour  Moe  surface 

STEF  ITER  OPTR  RAO  (INI 
371«  3  1  3.84T5 


63.A000  SECONDS  - 

HiMLt  -  HEDGE  -  HEAT  COEEE 

(BTU/LBI  (BTU/LBI  ILB/SQ  FT-SECt 

-IGOS.SB  -1436.27  l.$4S3 


- CHEMICAL  EROSION  ENERGY - 

- RATE  EBTU/SG  FT-SECI - 


0.0 


- NASS  LOSS  SUMMARY - 


- NASS  FtOR 

(LB/SQ  FT-SECl 

fCCCCSJltlM  fflllC 

(IN/ SEC  1 

RECESSION 

(INI 

T  MCftMOdICNl  vAL 

ABLATION  0.47500-01 

- ^  -  -  - _ 

0.47290-02 

0.10300F00 

W  1  E 

MECHANICAL 

EROSION  0.0 

0.0 

0.40410FOO 

CH/CHC 

0.96700 


I 
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